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EXECUTIVE  SUMMARY 


The  main  objective  of  this  work  is  to  develop  a  non-periodic  inspection  procedure  for  marine 
structures  so  as  to  maintain  their  reliability  at  a  prespecified  design  level  throughout  their  life. 
In  the  first  chapter  of  this  work,  a  summary  of  inspection  procedures  for  a  range  of  marine 
steel  structures  is  provided.  Specifically,  inspection  procedures  for  bulk  carriers,  fixed  offshore 
structures  and  semi-submersible  units  are  reviewed.  This  review  shows  that  there  is  very  little 
information  available  to  assess  the  reliability  of  these  inspection  techniques.  It  is  noted  that  only 
one  reliable  probability  of  detection  (POD)  curve  was  found.  In  view  of  the  limited  information 
available  to  assess  the  reliability  of  inspection  procedures  for  marine  steel  structures,  a  review 
of  POD  curves  used  in  the  aerospace  industry  is  presented  in  the  second  chapter.  This  review 
is  done  because  it  is  believed  that  the  shape  of  POD  curves  used  in  the  aerospace  industry  can 
provide  useful  guidelines  for  assessing  the  reliability  of  flaw  detection  and  at  the  same  time,  for 
establishing  POD  curves  for  marine  structures.  In  this  chapter,  emphasis  is  given  to  the  log  odds 
model  which  is  investigated  extensively.  Another  very  good  model  is  the  Weibull  model  which 
has  been  well  studied  in  dealing  with  applications  and  found  to  be  very  reasonable  for  the  POD 
function.  Indeed,  it  is  a  special  case  of  the  Weibull  model  that  is  used  as  a  POD  curve  in  the 
third  chapter.  In  the  third  and  main  chapter  of  this  report,  a  non-periodic  inspection  procedure 
is  developed  based  on  Bayesian  upgrading  and  taking  into  account  the  detailed  record  of  the 
entire  inspection  history  including  repair  or  replacement  records  for  each  and  every  component 
of  the  structure.  It  is  assumed  that  different  components  of  the  structure  are  subjected  to  different 
stress  levels.  In  the  research  done  in  the  past  using  Bayesian  analysis,  all  the  components  of  the 
structure  were  assumed  to  be  subjected  to  the  same  stress  level.  This  assumption  was  not  realistic 
and  therefore  different  stress  levels  are  considered  for  different  components  of  the  structure.  A 
numerical  example  is  provided  that  verifies  the  validity  and  effectiveness  of  Bayesian  analysis  to 
determine  appropriate  inspection  intervals  for  marine  structures  so  as  to  maintain  their  reliability 
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at  a  prespecified  design  level  throughout  their  life.  The  above-described  main  objective  of  this 
work  is  carried  out  in  the  thud  chapter. 

The  state-of-the-art  in  current  inspection  procedures  of  marine  structures  has  been  found  to 
be  generally  incompatible  with  the  probabilistic  approach  this  work  intends  to  develop.  Hence, 
the  first  two  chapters  describing  the  results  of  the  two  tasks  called  for  by  the  contract  could  not 
provide  probabilistic  information  directly  usable  in  Chapter  3.  This  fact  has  contributed  to  the 
appearance  that  the  first  three  chapters  of  this  report  are  somewhat  unrelated,  although  they  are 
all  an  integral  part  of  this  work. 

Finally,  further  study  is  suggested  in  the  following  five  areas: 

a.  In  this  work,  the  three  parameters  0,  c,  and  d  introducing  uncertainty  to  the  time  to  crack 
initiation,  fatigue  crack  propagation  and  probability  of  crack  detection  respectively,  were 
considered  to  be  uncorrelated.  However,  there  is  strong  evidence  that  /3,  c,  and  d  arc  in 
reality  correlated  to  each  other.  Therefore,  the  statistical  correlation  among  /?,  c,  and  d  and 
the  effect  of  this  correlation  on  the  obtained  results  require  further  study.  Another  aspect 
of  future  work  concerning  parameters  /?,  c,  and  d  is  to  examine  their  sensitivity  on  the 
obtained  results; 

b.  The  effect  of  the  form  of  certain  POD  curves  on  the  reliability  of  marine  structures  throughout 
their  service  life  requires  further  study.  A  comparison  has  to  be  made  among  several 
established  POD  curves  in  the  aerospace  industry,  in  order  to  assess  their  relative  influence 
on  the  reliability  of  marine  structures  subjected  to  non-periodic  inspections.  In  this  way, 
more  reliable  POD  curves  can  be  established  for  marine  structures; 

c.  The  cost-effectiveness  of  the  proposed  method  of  non-periodic  inspections  based  on  Bayesian 
analysis  requires  further  study.  Specifically,  a  cost-benefit  analysis  can  be  performed  taking 
into  consideration  the  cost  of  the  non-periodic  inspection  procedure  and  the  increased  level 
of  reliability  for  the  structure.  These  results  have  to  be  compared  with  the  results  of  the 
cost-benefit  analysis  associated  with  the  standard  periodic  inspection  procedure; 

d.  The  verification  of  the  proposed  methodology  using  actual  data  from  inspections  of  marine 
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structures  is  one  of  the  most  important  tasks  of  future  work.  This  task  can  be  accomplished  by 
taking  advantage  of  already  completed  inspections  of  marine  structures  to  determine  whether 
these  structures  actually  maintained  a  prespecified  reliability  level  throughout  their  life; 
e.  The  failure  rate  expression  after  crack  initiation  should  at  least  be  validated  by  Monte  Carlo 
simulation  utilizing  die  crack  propagation  law  and  uncontrolled  crack  growth  condition  based 
cm  fracture  mechanics  theory  under  various  random  stress  histories  consistent  with  the  stress 
intensity  factor  fluctuation  considered. 
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SUMMARY 


A  summary  of  inspection  procedures  for  a  range  of  marine  steel  structures  is  provided  in 
the  first  chapter  of  this  work.  Specifically,  inspection  procedures  for  bulk  carriers,  fixed  offshore 
structures  and  semi-submersible  units  are  reviewed.  This  review  shows  that  there  is  very  little 
information  available  to  assess  the  reliability  of  these  inspection  techniques.  It  is  noted  that  only 
one  reliable  probability  of  detection  (POD)  curve  was  found.  In  view  of  the  limited  information 
available  to  assess  the  reliability  of  inspection  procedures  for  marine  steel  structures,  a  review 
of  POD  curves  used  in  the  aerospace  industry  is  presented  in  the  second  chapter.  This  review 
is  done  because  it  is  believed  that  the  shape  of  POD  curves  used  in  the  aerospace  industry  can 
provide  useful  guidelines  for  assessing  the  reliability  of  flaw  detection  and  at  the  same  time, 
for  establishing  POD  curves  for  marine  structures.  In  this  chapter,  emphasis  is  given  to  the  log 
odds  model  which  is  investigated  extensively.  Another  very  good  model  is  the  Weibull  model 
which  has  been  well  studied  in  dealing  with  applications  and  found  to  be  very  reasonable  for 
the  POD  function.  Indeed,  it  is  a  special  case  of  the  Weibull  model  that  is  used  as  a  POD 
curve  in  the  third  chapter. 

The  third  chapter  constitutes  the  main  part  of  this  work  whose  basic  objective  is  to  develop 
a  non-periodic  inspection  procedure  for  marine  structures  so  as  to  maintain  their  reliability  at  a 
prespecified  design  level  throughout  their  life.  This  procedure  is  based  on  Bayesian  upgrading 
and  takes  into  account  the  detailed  record  of  the  entire  inspection  history  including  repair  or 
replacement  records  for  each  and  every  component  of  the  structure.  It  is  considered  that 
different  components  of  the  structure  are  subjected  to  different  stress  levels.  In  the  research 
done  in  tire  past  using  Bayesian  analysis,  all  the  components  of  the  structure  were  assumed  to 
be  subjected  to  the  same  stress  level.  This  assumption  was  not  realistic  and  therefore  different 
stress  levels  are  considered  for  different  components  of  the  structure.  Finally,  a  numerical 
example  is  provided  that  verifies  the  validity  and  effectiveness  of  Bayesian  analysis  to  determine 
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appropriate  inspection  intervals  for  marine  structures  so  as  to  maintain  their  reliability  at  a 
prespecified  design  level  throughout  their  life. 
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I.  REVIEW  OF  CURRENT  INSPECTION  PROCEDURES 
OF  MARINE  STRUCTURES 


t.t  Introduction 

This  chapter  presents  results  of  work  performed  by  Stewart  Technology  Associates  (STA), 
subcontractor  to  Modem  Analysis,  Inc.  on  a  project  for  the  US  Coast  Guard. 

This  project  deals  with  the  reliability  of  marine  structures,  concentrating  on  fatigue  damage 
and  its  detection  during  regular  in-service  inspection.  The  main  thrust  of  the  work  performed 
by  STA  was  to  summarize  inspection  procedures  for  a  range  of  marine  steel  structures  and  to 
provide  summaries  of  inspection  findings,  in  cooperation  with  classification  societies.  The  work 
was  directed  towards  the  structural  integrity  of  the  main  hull,  or  main  structure,  of  each  of  the 
marine  structure  categories  considered.  Fatigue  damage,  as  evidenced  by  surface  cracks,  was 
the  principal  type  of  damage  to  be  considered  resulting  in  a  reduction  in  strengths. 

STA  visited  the  American  Bureau  of  Shipping  (then)  in  New  York,  and  Det  norske  Veritas 
in  Norway,  in  order  to  discuss  the  inspection  procedures  and  results.  Additionally,  three  visits 
were  made  to  key  individuals  and  organizations  in  London,  as  well  as  to  Exxon  in  Houston. 
Telephone  discussions  with  other  companies  also  contributed  to  the  general  picture  of  industry 
experience  and  current  practice  presented  in  these  earlier  reports. 

This  report  summarizes  the  key  information  presented  in  earlier  reports  and  additionally 
presents  some  further  key  published  information  relating  to  the  probability  of  detecting  cracks 
underwater  when  marine  structures  are  inspected.  During  the  course  of  the  work  by  STA,  it 
became  clear  that  the  industry  had  very  little  information  enabling  any  kind  of  assessment  of 
the  reliability  of  inspection  techniques  to  be  made.  Consequently,  STA  was  directed  by  Modem 
Analysis  to  try  to  obtain  any  available  probability  of  detection  (POD)  information.  Only  one 
reliable  POD  curve  was  found,  and  this  was  in  cooperation  with  DnV. 
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1.2  Summary  of  Inspection  Procedures 

1.2.1  Bulk  Carriers 

Both  the  American  Bureau  of  Shipping  (ABS)  and  Det  norske  Vferitas  (DnV)  have  rule* 
that  govern  the  inspection  requirements  for  bulk  carriers,  including  dry  bulk  and  oil  tankers. 
These  requirements  concentrate  principally  upon  corrosion  and  thickness  testing.  All  inspection 
for  cracks  is  visual  unless  there  is  something  unusual,  and  then  it  is  up  to  the  judgment  of  the 
individual  surveyor.  For  tankers,  particularly  in  the  DnV  rules,  the  inspection  requirements  are 
greater  than  for  dry  cargo  vessels,  but  there  is  still  no  requirement  for  anything  more  than  visual 
inspection  to  locate  cracks.  Thickness  measurements  are  required  at  special  periodic  surveys. 
These  special  periodic  surveys  are  typically  at  2-1/2  or  4-year  intervals.  No  guidance  is  given 
in  these  rules  as  to  acceptable  defect  sizes  if  cracks  are  found.  It  is  up  to  the  judgment  of  the 
individual  surveyor  as  to  whether  or  not  the  cracks  will  be  repaired. 

All  inspection  normally  takes  place  with  ships  in  the  dry.  This  is  done  by  internal  inspection 
of  the  ship’s  structure,  either  while  die  ship  is  floating  or  dry -docked.  External  inspection  of  the 
hull  structure  is  normally  done  in  dry  dock.  Both  classification  societies  have  some  provision 
for  underwater  surveys,  but  both  eventually  require  dry-docking  after  a  number  of  years. 

Neither  classification  society  has  adopted  a  philosophy  for  calculating  the  growth  rate  of 
defects  that  may  develop  in  ship  hull  structures.  This  is  principally  because  through-hull  cracks 
are  generally  detected  because  of  leakage  into  cargo  or  other  normally  dry  spaces,  and  in  the 
event  of  a  through-hull  crack,  it  is  repaired  immediately. 

1.2.2  Fixed  Offshore  Structures 

In  1986  a  notice  of  proposed  rule-making  was  published  in  the  Federal  Register.  This 
proposed  an  inspection  requirement  for  structures  on  the  U.S.  offshore  continental  shelf  (OCS). 
The  requirement  was  rather  general,  requiring  that  periodic  inspection  of  such  structures  is 
performed  to  determine  “the  condition  of  the  entire  structure.”  An  annual  report  from  die  operator 
is  to  be  submitted  to  the  MMS,  “stating  which  platforms  have  been  inspected  in  the  preceding 
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12  months,  the  extent  and  area  of  inspection,  and  the  type  of  inspection  employed,  i.e.,  visual, 
magnetic  particle,  ultrasonic  testing.”  In  the  absence  of  existing  requirements  to  inspect  OCS 
structures,  it  is  noted  that  operators  perform  their  own  inspections,  sometimes  within  a  carefully 
planned  framework  of  life-cycle  costing,  but  more  frequently  without  an  overall  plan.  Inspection 
procedures  vary  from  simply  checking  periodically  that  the  structure  is  still  there  (unmanned 
structures)  to  cleaning  and  NDT  of  critical  joints,  underwater,  on  a  regular  planned  basis.  It  is 
true  to  say  that  the  level  of  inspection  of  fixed  structures  in  OCS  waters  is  considerably  lower 
than  that  undertaken  for  structures  in  the  North  Sea. 

For  structures  in  the  North  Sea,  governmental  requirements  are  more  detailed,  and  inspection 
procedures  are  more  rigorous  and  certainly  cost  a  great  deal  of  money.  Every  year  a  considerable 
amount  of  in-service  inspection  of  offshore  structures  in  North  Sea  waters  is  performed  in  order 
to  ensure  the  safety  of  personnel  and  production.  The  Norwegians,  in  particular,  have  developed 
systematic  methods  based  on  probabilistic  models  and  cost  resource  allocation  for  the  inspection 
of  structures  in  Norwegian  waters. 

DnV  have  produced  rules  addressing  the  questions  of  personnel  qualifications,  inspection 
procedures,  and  equipment  capability.  These  rules  also  differentiate  between  three  types  of 
inspection: 

Type  I  (Green): 

General  visual  inspection  to  detect  obvious  damage.  Prior  cleaning  of  inspection 
items  is  not  needed. 

Type  II  (Blue): 

Close  visual  inspections  to  detect  hidden  damage.  Prior  cleaning  of  inspection  items 
is  normally  necessary. 

Type  ffl  (Red): 
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Close  visual  inspection  and  testing  (NDT  or  DT)  to  detect  incipient  or  hidden 
damage.  Prior  cleaning  of  inspection  items  is  required. 

DnV  cite  two  basic  methods  for  planning  an  in-service  inspection  program,  one  based  upon 
the  design,  fabrication,  and  installation  (DFI)  knowledge,  and  the  other  based  upon  the  struc¬ 
ture’s  condition  records  (SCR).  DFI  principally  identifies  areas  where  inspection  should  be  con¬ 
centrated,  based  upon  calculation  of  the  structure's  predicted  performance  in  the  ocean  environ¬ 
ment,  taking  account  of  any  deviations  from  the  original  design  occurring  during  fabrication  and 
installation.  SCR  is  a  system  of  continuous  revision  to  take  account  of  inspection  findings  in 
service,  when  it  may  be  found  in  practice  that  deterioration  of  some  areas  occurs  more  rapidly 
than  predicted  with  DFI  as  the  inspection  basis.  DnV  emphasize  the  bookkeeping  aspects  of 
inspection  results,  enabling  trend  analyses  to  be  relatively  easily  undertaken.  As  with  the  ABS, 
they  refer  to  an  inspection  program,  specific  to  each  installation;  but  the  DnV  rules  give  much 
guidance  as  to  the  form  of  such  a  program,  while  the  ABS  rules  do  not.  DnV  rules  also  offer 
guidance  as  to  selection  or  areas,  in  general,  for  inspection  and  offer  possible  types  and  causes 
of  defects  that  may  be  found. 

One  of  the  most  useful  pieces  of  information  giving  a  picture  of  inspection  findings  for  fixed 
steel  structures  in  the  Norwegian  Sector  of  the  North  Sea  was  a  table  provided  by  DnV  for  piled 
steel  structures  they  have  inspected  in  the  years  1975  to  1984.  This  table  is  reproduced  below; 

Table  1.1  PILED  STEEL  STRUCTURES 
1975-1984 

TOTAL  NOS:  21 


Nos.  of  Platforms  With  Defects 

WIRE  SCAFING 

18 

DEFLECTED  MEMBERS 

11 

DENTED  MEMBERS 

15 

9 


Table  1.1  (Continued)  PILED  STEEL  STRUCTURES 

1975-1984 


TOTAL  NOS:  21 


MISSING,  LOOSE,  DEFECT  ANODES 

14 

GENERAL  CORROSION 

2 

PITTING  CORROSION 

20 

BURN  MARK 

4 

HEAVY  MARINE  GROWTH 

13 

SCOUR 

1 

DEBRIS 

21 

CONFIRMED  CRACK 

12 

PROPAGATING  CRACK 

3 

As  can  be  seen  from  this  table,  in  all  these  years  of  inspection  on  the  21  structures,  only 
three  propagating  cracks  were  found  (in  major  structural  members,  not  secondary  members  such 
as  conductor  bay  framing). 

Figure  1.1  seems  to  represent  the  best  POD  information  available  in  the  industry  at  1987, 
which  indeed  will  be  used  in  Chapter  II  where  a  method  is  developed  to  more  rigorously 
estimate  POD  functions  on  a  similar  database. 
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Figure  1.1  Example  of  Crack  Detection  Probabilities  for  Magnetic  Particle  Testing  Under 
Water  Based  on  Measurements  by  Moncaster.  Included  Are  200  Observations  on  14  Cracks 

Although  very  few  propagating  cracks  had  been  found  at  that  time  on  main  structural 
members,  there  was  a  lot  of  information  available  about  cracks  found  on  conductor  bay  framing, 
particularly  as  a  result  of  design  errors  in  the  1970’s.  These  design  errors  had  omitted  considering 
vertical  wave  force  effects  on  the  conductor  bay  framing.  These  vertical  forces  resulted  in  bounce 
action  of  the  bay  framing  and  often  led  to  rapid  fatigue  failures.  Although  potentially  dangerous 
in  that  risers  would  lose  support,  the  overall  structural  integrity  was  not  affected  by  loss  of 
conductor  bays. 

1.23  Semi-Submersible  Units 

The  ABS  requirements  for  inspection  of  semi-submersibles  in  ABS  class  are  given  in  the  ABS 
Rules  for  Building  and  Classing  Mobile  Offshore  Drilling  Units.  The  minimum  requirements,  as 
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stated  in  the  rules,  are  for  annual  and  4-yearly  inspections  (special  periodic  surveys).  The  rules 
also  provide  for  continuous  surveys,  and  this  is  becoming  the  most  common  practice  in  industry. 
Dry-docking  is  required  every  2-1/2  years.  Recently  the  ABS  have  emphasized  the  advantage 
of  having  an  inspection  manual  for  each  rig.  Reference  1.1  contains  their  suggestions  for  an 
inspection  manual  as  presented  at  a  February  1989  meeting  of  the  Society  of  Naval  Architects 
and  Marine  Engineers  in  Houston. 

As  noted  in  die  original  STA  report,  the  growth  or  fatigue  cracks  in  semi-submersible  units 
has  been  the  subject  of  much  research  in  the  offshore  industry.  Units  operating  in  the  Gulf  of 
Mexico  and  other  relatively  calm  water  areas  do  not  suffer  the  fatigue  damage  to  which  units 
operating  in  harsh  environments  are  prone.  Since  the  1980  Alexander  L.  Kielland  sinking  in  a 
storm  in  the  North  Sea,  when  123  lives  were  lost,  a  great  deal  of  attention  has  been  centered  upon 
detecting  fatigue  cracks  and  predicting  the  rate  of  growth  of  fatigue  cracks  in  harsh  environment 
semi-submersibles.  DnV  have  evolved  inspection  procedures  for  semi’s  involving  a  monthly 
inspection  of  all  critical  joints  for  through-thickness  cracks  that  can  be  internally  detected  by 
leakage  into  normally  dry  spaces.  Their  philosophy  is  that  a  close  visual  inspection  of  all  critical 
joints  will  be  undertaken,  using  MPI,  every  four  to  five  years.  However,  they  believe  the  time 
taken  for  a  crack  to  propagate  to  through-thickness  is  around  two  to  five  years  in  the  North  Sea 
for  typical  semi-submersibles.  Hence  they  have  an  intermediate  survey  every  two  and  one-half 
years,  or  they  have  monthly  inspection  for  leaks  into  normally  dry  areas.  They  estimate  that 
the  time  for  a  through-thickness  crack  to  propagate  until  member  separation  is  typically  several 
months.  Additionally,  they  have  now  introduced  a  redundancy  requirement  such  that  rigs  must 
be  able  to  survive  at  least  a  one-year  storm  after  the  failure  of  any  individual  bracing  or  bracing 
connection. 

It  was  clear  from  both  DnV’s  experience  with  extraordinary  surveys  of  rigs  following 
the  Alexander  Kielland  accident  and  from  discussions  with  Exxon,  who  had  been  performing 
their  own  special  surveys  of  rigs,  that  many  cracks  were  missed  during  in-service  inspections. 
Furthermore,  these  cracks  had  often  been  missed  during  several  in-service  inspections  and  were. 
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in  fact,  attributed  to  fabrication  defects  which  had  existed  since  the  rig  left  the  fabrication  yards. 
Although  this  knowledge  was  common  in  industry,  there  was  almost  no  quantification  of  the 
probability  of  cracks  being  missed  during  in-service  inspections.  It  should  also  be  noted  that 
these  cracks  were  missed  on  rigs  which  were  being  inspected  in  relatively  good  conditions,  i.e., 
benign  weather,  above  the  water,  without  the  need  for  special  rigging. 

1.3  New  Information 

Reference  1.2  is  a  paper  presented  at  OMAE  Europe  89’  in  March  1989  by  Inge  Lotsberg  and 
Finn  Kiricemo  of  DnV.  This  paper  describes  a  systematic  method  of  planning  in-service  inspection 
of  steel  offshore  structures.  The  method  centers  upon  a  mathematical  model  for  optimization  of 
in-service  inspection  based  upon  reliability  methods  and  resource  allocation.  Bayesian  updating 
of  reliability  is  included  in  the  model.  This  paper  includes  the  same  probability  of  detection 
(POD)  curve  that  was  produced  in  the  STA  initial  repents.  It  is  assumed  in  the  paper  that 
inspection  for  fatigue  cracks  below  water  is  performed  using  magnetic  particle  inspection  (MPI). 
It  is  a  POD  curve  for  MPI  that  is  produced  in  the  paper.  It  is  also  noted  in  the  paper  that  the 
probability  of  detecting  a  crack  is  a  function  of  its  surface  length.  Accounting  for  uncertainty  in 
fatigue  life  predictions  and  reliability  of  the  inspection  method  in  the  mathematical  model,  the 
fatigue  reliability  in  the  paper  is  updated  according  to  Bayes  rule.  Inspection  intervals  can  then 
be  optimized  while  maintaining  a  prescribed  minimum  reliability  level.  Due  to  this  updating  of 
reliability,  the  inspection  interval  is  increased  as  a  function  of  service  life.  This  result  is  based 
upon  no  crack  being  detected  during  inspection. 

Another  conclusion  from  the  paper  is  that  it  is  cost-efficient  to  base  die  inspection  reliability 
on  a  crack  depth  equal  to  one-third  of  the  thickness  of  the  joint,  due  to  a  low  repair  cost  for  this 
crack  depth  compared  to  a  crack  depth  completely  through  the  thickness  of  the  joint  This  is 
because  a  crack  depth  of  just  one-third  of  the  joint  thickness  can  be  repaired  by  grir  ding.  The 
repair  for  a  crack  that  has  gone  through  the  thickness  of  the  joint  requires  underwater  welding 
and  is  therefore  enormously  more  expensive. 
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It  is  also  noted  that  the  authors  of  this  paper  assumed  that  an  aspect  ratio  for  crack  depth 
to  length  of  0.1S  is  appropriate  to  the  cracks  and  joints  of  offshore  structures  that  the  paper 
addresses.  Their  POD  curves  gives  a  90-percent  probability  of  detecting  a  3-millimeter  deep 
crack. 

Reference  1.3  is  another  paper  presented  at  the  March  OMAE  conference.  This  paper  is 
by  Paul  Frieze  and  Jacob  Kam.  Its  title  is  ‘The  Reliability  Assessment  of  the  Nondestructive 
Inspection  of  Offshore  Tubular  Structures.”  In  this  paper  the  authors  note  the  historical  aerospace 
inspection  requirement  to  achieve  a  90/95  success  rate.  That  is,  inspection  should  achieve 
90-percent  POD  at  the  95-percent  level  of  confidence.  These  figures  are  derived  assuming 
POD  trial  results  can  be  statistically  quantified  by  nomial  distribution.  This  approach  has  been 
basically  inherited  by  the  offshore  industry.  However,  in  order  to  achieve  this  level  of  success, 
28  cracks  in  a  particular  range  are  all  required  to  be  detected.  If  one  is  missed,  46  successful 
detections  are  necessary  if  the  same  level  of  success  is  to  be  realized.  It  has  been  noted  recently 
that  should  success  fall  significantly  below  this  level,  the  relevant  statistical  distribution  may  no 
longer  be  binomial  but  some  unknown  distribution. 

The  authors  quantify  POD  success  using  a  variety  of  geometrical  definitions.  These  are 
as  follows: 

•  Size  ratio  (or  length  ratio) 

•  Coverage  ratio 

•  Overlap  ratio 

Size  (length)  ratio  is  defined  as  the  ratio  4^,  where  LM  and  LA  are  the  measured  and  the 
actual  crack  sizes  (lengths),  respectively.  Values  can  range  from  zero  to  infinity.  Zero  indicates 
a  miss,  infinity  a  spurious  indication,  and  unity  a  completely  successful  indication.  With  an 
average  inspection  technique,  plots  of  frequency  versus  length  ratio  are  shown  in  Fig.  1.2. 
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Fig.  1.2  Typical  distributions  of  inspection  results  based  on  length 
ratio. 


POD  curves  can  then  be  generated  for  any  “threshold”  value  of  length  ratio.  These  POD 
curves  represent  the  area  under  the  curves  in  Figure  1.2  to  the  right  of  the  threshold  level.  Figure 
1.3  illustrates  POD  curves  for  0,  50  and  90  percent  threshold  levels.  They  all  converge  to  0 
POD  as  crack  length  diminishes,  as  few  techniques  are  successful  in  this  range. 


Actual  length 


Fig.  1.3  Typical  length  ratio  POD  curves. 


Unfortunately,  the  authors  are  constrained  by  confidentiality  agreements  with  their  clients  and 
have  not  put  numbers  to  the  actual  length  axis  of  their  POD  curve.  However,  it  should  be  noted 
that  the  general  shape  of  the  curve  is  similar  to  that  produced  by  STA  in  their  earlier  reports. 

The  reader  of  this  report  is  encouraged  to  refer  to  Ref.  1.3  to  understand  the  authors’ 
definitions  of  coverage  ratio  and  overlap  ratio.  In  the  view  of  STA,  both  these  definitions  are 
helpful  in  understanding  the  significant  differences  between  attempts  to  define  POD  curves  for 
the  offshore  industry  and  POD  curves  that  have  historically  been  established  in  the  aerospace 
industry. 

Another  important  aspect  of  Ref.  1.3  is  the  cost  information  that  it  presents  relating  to  un¬ 
derwater  inspection  and  repair.  A  diving  team  often  consists  of  eight  to  twelve  divers/inspectors, 
which  for  nonsaturation  diving  currently  costs  some  $7,000  per  day.  For  a  short  run  of  two  to 
three  weeks,  the  total  cost  is  nearly  $140,000.  Usually  two  shifts  are  used  to  make  maximum 
use  of  the  weather  window.  Full  saturation  diving  is  considerably  more  expensive.  A  diving 
support  vessel  costs  around  $50,000  per  day.  With  mobilization  costs,  specialist  equipment  hire, 
and  daily  supplies,  a  modest  full  saturation  inspection  generates  a  minimum  outlay  in  the  region 
of  $1.4  million.  For  year-round  surveys  of  a  large  offshore  field  in  the  North  Sea,  an  operator 
may  have  to  spend  as  much  as  $24  million.  The  authors  point  out  that  there  is  presently  no 
objective  criterion  to  determine  how  successful  the  surveys  are  or  what  they  have  achieved  for 
this  enormous  investment. 
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References  1.4  and  1.5  are  two  articles  from  Ocean  Industry,  March  1987.  The  first  one 
by  W.P.  Stewart  [1.4]  includes  the  original  POD  curve  presented  in  the  initial  STA  report  and 
reproduced  here  in  Fig.  1.4. 


Fig.  1.4  Shaded  area  shows  results  from  a  UK  DOE/Shell  North  Sea  study  where  MPI  was 
applied  underwater  to  detect  fatigue  cracks.  A  95°/o  confidence  band  is  shown  for  the  probabil- 
Ity-of-detection  regression  curve.  A  total  of  34  defects  was  inspected,  with  342  observations. 
The  colored  curve  in  the  plot  shows  results  from  in-air  inspection  of  aircraft  structures;  the 
mean  value  curve  is  shown  from  a  series  of  60  inspections  of  more  than  20  defects,  with  over 
2,000  observations.  The  lower  black  curve  is  a  more  conservative  viewpoint  of  typical  POD 
values  commonly  experienced  for  above-water  inspection  of  rigs  when  MPI  is  used. 


The  second  [1.5]  is  a  paper  by  F.R.  Frisbie  on  “Inspecting  and  Repairing  Offshore  Platforms 
Today.”  This  paper  gives  some  interesting  cost  information  in  Figure  1.5. 
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from  the  U.K.  Department  of  Energy  published  in  1984.  It  is  included  since  it  was  sent  to 
Dr.  Shinozuka  at  Modem  Analysis  by  Mr.  Ken  Bitting  from  the  U.S.  Coast  Guard  R&D 
Center  following  a  request  from  the  SR- 13 17  Project  Technical  Committee  members  for  any 
information  on  probability  of  detection  data  for  underwater  inspection  systems.  The  information 
this  document  contains  is  of  some  interest  in  that  it  includes  data  using  ultrasonic  measurements 
and  MPI  measurements  underwater.  However,  it  principally  shows  some  of  the  difficulties  with 
correlating  this  type  of  data  and  the  almost  impossible  task  of  using  this  type  of  data  to  develop 
POD  curves. 

1.4  Conclusions 

In  this  chapter,  a  summary  of  inspection  procedures  for  a  range  of  marine  steel  structures 
was  provided.  Specifically,  inspection  procedures  for  bulk  carriers,  fixed  offshore  structures  and 
semi-submersible  units  were  reviewed.  This  review  showed  that  there  is  very  little  information 
available  to  assess  the  reliability  of  these  inspection  techniques.  It  is  noted  that  only  one  reliable 
probability  of  detection  (POD)  curve  was  found.  In  view  of  this  limited  information  available 
to  assess  the  reliability  of  inspection  procedures  for  marine  steel  structures,  a  review  of  POD 
curves  used  in  the  aerospace  industry  will  be  presented  in  the  following  chapter.  This  review 
will  be  done  because  it  is  believed  that  the  shape  of  POD  curves  used  in  the  aerospace  industry 
can  provide  useful  guidelines  for  assessing  the  reliability  of  flaw  detection  and  for  establishing 
POD  curves  for  marine  structures. 
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.  RELIABILITY  OF  FLAW  DETECTION 

2.1  Introduction 

In  view  of  the  limited  information  available  to  assess  the  reliability  of  inspection  procedures 
for  marine  steel  structures  as  concluded  in  the  first  chapter,  a  review  of  POD  curves  used  in 
the  aerospace  industry  is  presented  in  this  chapter.  This  review  is  done  because  it  is  believed 
that  the  shape  of  POD  curves  used  in  the  aerospace  industry  can  provide  useful  guidelines  for 
assessing  the  reliability  of  flaw  detection  and  for  establishing  POD  curves  for  marine  structures. 

Current  nondestructive  inspection  (NDI)  systems  are  not  capable  of  repeatedly  producing 
correct  indications  when  applied  to  flaws  of  the  same  length.  The  chance  of  detecting  a  given 
crack  length  depends  on  many  factors,  such  as  the  location,  orientation  and  shape  of  the  flaw, 
materials,  inspectors,  inspection  environments,  etc.  As  a  result,  the  probability  of  detection 
(POD)  for  all  cracks  of  a  given  length  has  been  used  in  the  literature  to  define  the  capability  of 
a  particular  NDI  system  in  a  given  environment. 

In  aerospace  applications,  a  nondestructive  inspection  limit,  a^DE,  is  chosen,  which  is 
a  crack  length  that  usually  corresponds  to  a  high  detection  probability  and  high  confidence 
level.  For  instance,  the  damage  tolerance  specifications  for  aircraft  structures  require  that  the 
NDI  system  be  capable  of  detecting  a  specified  crack  length,  qnde ,  at  a  particular  location 
with  a  90%  detection  probability  and  95%  confidence  level  (see  Ref.  2.1).  The  fracture 
mechanics  propagation  life,  Nf,  is  the  life  for  crack  length  o^de  to  propagate  to  the  critical 
crack  length  ac,  under  expected  usage  environments.  The  inspection  interval,  r,  is  equal  to 
Nf  divided  by  a  safety  factor,  S/,  i.e.,  t  =  N//S/.  In  evaluating  the  structural  reliability  under 
scheduled  inspection  maintenance,  however,  the  information  of  cl^de  is  of  little  value  and  the 
uncertainty  of  the  NDI  system  should  be  taken  into  account  [2.2— 2.5]. 

Flaw  detection  reliability  is  defined  as  the  probability  of  detecting  a  flaw  under  pre-specified 
inspection  conditions.  This  probability  is  a  function  of  the  crack  length.  Figure  2.1  shows  a  plot 
of  inspection  results  for  individual  cracks  emanating  from  fastener  holes  in  a  skin  and  stringer 
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wing  assembly  using  eddy  current  surface  scans  [2. 6-2.7].  The  points  represent  the  proportion 
of  times  individual  flaws  were  detected  versus  the  length  of  the  flaw.  This  figure  illustrates  that 
although  the  detection  probabilities  of  individual  flaws  generally  increase  with  crack  length,  not 
all  flaws  of  the  same  length  have  the  same  detection  probability.  This  variability  in  detection 
probabilities  at  a  crack  length  requires  a  consistent  definition  of  the  probability  of  detection 
(POD)  as  a  function  of  crack  length  “a.”  The  POD  (a)  function  is  defined  as  the  proportion 
of  flaws  that  will  be  detected  as  a  function  of  crack  length,  i.e.,  the  mean  trend  in  detection 
probabilities  as  a  function  of  crack  length  [2.6-2.7].  The  solid  curve  in  Fig.  2.1  is  a  POD 
function  obtained  from  the  inspection  results  (points)  as  will  be  described  later. 

It  follows  from  Fig.  2.1,  that  an  NDI  system  may  result  in  two  types  of  incorrect  indications 
(i)  failure  to  give  a  positive  indication  in  the  presence  of  a  crack  whose  length  is  greater  than 
aNDE*  referred  to  as  a  Type  I  error,  and  (ii)  give  a  positive  indication  when  the  crack  length  is 
smaller  than  a^DE^  referred  to  as  a  Type  II  error.  The  Type  I  error  allows  components  containing 
a  crack  length  longer  than  a^oE  to  remain  in  service,  thus  greatly  increasing  the  potential  safety 
hazard.  For  safety  critical  components,  the  Type  I  error  is  of  primary  concern.  The  Type  II  error 
rejects  good  components  and,  hence  has  an  adverse  effect  on  the  cost  of  repair/replacement  and 
life  cycle  cost.  In  applications  such  as  retirement-for-cause  (RFC)  life  management,  however, 
both  Type  I  and  Type  II  errors  are  important,  because  the  criterion  used  in  RFC  life  management 
is  the  minimization  of  the  life-cycle-cost  (LCC)  [2.8-2.9J.  For  a  given  NDI  system  with  a  single 
inspection,  it  is  impossible  to  reduce  the  Type  II  error  without  increasing  the  Type  I  error  and 
vice  versa  [2.10-2.12].  It  is  obvious  that  the  ideal  inspection  capability  of  an  NDI  system  is  a 
unit  step  function.  Figure  2.2  schematically  shows  an  ideal  and  a  realistic  POD  function.  The 
ideal  inspection  system  would  detect  all  flaws  larger  than  a^DE  and  none  smaller  than  a^DE 
as  indicated  by  a  unit  step  function  in  Fig.  2.2,  in  which  both  Type  I  and  Type  II  errors  are 
zero.  Unfortunately,  such  an  ideal  NDI  system  is  far  from  reality.  Technical  approaches  to 
reduce  both  types  of  errors  using  multiple  inspection  procedures  were  studied  recently  by  Yang 
et  al.  [2.10-2.12]. 
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Many  factors  influence  the  capability  of  an  NDI  system  to  identify  flaws  in  a  structure. 
These  include  (1)  system  factors  which  affect  the  ability  of  the  system  to  consistently  produce 
and  interpret  the  information  upon  which  flaw  decisions  are  made,  and  (2)  factors  which  are 
characteristics  of  the  individual  flaws  being  inspected.  The  net  effect  is  uncertainty  in  the 
detection  of  flaws  so  that  the  process  of  quantifying  the  inspection  capability  of  a  particular 
system  requires  a  careful  NDI  reliability  demonstration  program  coupled  with  a  probabilistically 
based  analysis  of  the  data.  This  section  describes  the  reliability  of  flaw  detection  and  the  analysis 
of  NDI  reliability  data  following  Refs.  2.6  and  2.7. 

As  described  previously,  the  detection  probability  for  a  given  crack  length  involves  consid¬ 
erable  statistical  variability.  The  distribution  of  detection  probabilities  at  a  given  crack  length 
is  illustrated  in  Fig.  2.3.  The  curve  connecting  the  average  values  of  the  detection  probabili¬ 
ties  for  all  crack  lengths  is  defined  as  the  POD  (a)  function.  Hence,  the  POD  (a)  function  is 
a  function  which  passes  through  the  mean  of  detection  probabilities  at  each  crack  length,  i.e., 
a  regression  function.  Consequently,  many  individual  cracks  will  have  detection  probabilities 
below  the  POD  (a)  value. 

2.2  POD  Functional  Form 

The  information  on  POD  functions  for  NDI  systems  is  needed  in  the  reliability  analysis  of 
structures  under  scheduled  inspection  maintenance.  It  is  also  crucial  for  the  determination  of  the 
inspection  interval.  To  establish  the  POD  function  from  experimental  test  results,  a  functional 
form  should  be  assumed.  The  so-called  log  odds  or  log  logistic  model  has  been  investigated 
extensively  [2.6-2.7], 

POD(«)  =  {l+«p[-|('^)]}"  (2.U) 

in  which  POD  (a)  is  the  probability  of  detecting  crack  size  a,  and  n  and  a  are  parameters. 
Methods  for  estimating  the  parameters  n  and  a  from  NDI  reliability  data  is  a  major  subject  of 
this  section. 
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Let  ao.s  be  the  median  crack  detection  capability,  i.e.,  the  crack  length  associated  with  a 
50%  detection  probability,  POD  (ao.s)  =  0.5.  Then,  it  follows  from  Eq.  2.1a  that 


H  =  lnao.5 


(2.1b) 


Thus,  the  parameter  pi  represents  the  central  location  of  the  POD  curve.  The  parameter  a  is 
a  measure  of  the  flatness  of  the  POD  function,  the  larger  the  value  of  <7,  the  slower  the  POD 
function  approaches  one.  The  parameters  pi  and  a  are  referred  to  as  the  location  and  scale 
parameters,  respectively. 

Another  POD  function,  referred  to  as  the  Weibull  function,  has  also  been  used  [2.2,  2.3, 

2.10]: 

POD  (a)  =  0  a  <  c 


in  which  e  is  the  crack  length  below  which  a  crack  cannot  be  detected  by  the  NDI  system. 
Again,  a  and  0  are  constants,  representing  the  bandwidth  and  central  location  of  the  POD 
function,  respectively. 

As  mentioned  previously,  the  POD  function  is  a  unit  step  function  at  asDE  for  an  ideal 
NDI  system,  i.e.. 


POD  (a)  ==  0 


a  <  onde 


(2.3) 


=  1  a  >  affDE 

Such  an  ideal  POD  function  can  be  obtained  from  the  Weibull  function  by  setting  e  = 
a  pipe,  0  -*  0  and  a  — ►  oo.  Hence  the  Weibull  function  given  by  Eq.  2.2  includes  the  ideal 
POD  curve,  Eq.  2.3,  as  a  special  case. 

These  and  other  POD  functions  have  been  proposed  in  the  literature  [e.g.,  2.2~2.5,  2.13, 
2.14].  Among  these,  the  Weibull  model  and  log  odds  model  appear  to  be  most  viable  for  die 
analysis  that  is  to  be  performed  here.  However,  the  Weibull  model  has  been  well  studied  in 
dealing  with  other  applications.  Therefore,  only  the  log  odds  model  has  been  investigated  here 
extensively  using  available  NDI  reliability  data  [2. 6,2.7]  and  it  has  been  shown  that  the  log  odds 
model  is  also  very  reasonable  for  the  POD  function.  In  Chapter  III,  however,  a  special  case  of 
the  Weibull  model  with  a  =  1  and  e  =  0  will  be  used. 
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2.3  Statistical  Estimation  of  POD  Function 


In  order  to  establish  the  POD  function  for  a  particular  NDI  system  associated  with  a  particular 
structural  detail  under  a  particular  inspection  environment,  NDI  reliability  demonstrations  should 
be  conducted.  Data  collected  from  NDI  reliability  demonstration  programs  consist  of  two 
categories:  (i)  data  in  which  the  inspection  result  (pass  or  fail)  is  recorded,  and  (ii)  data  in 
which  the  response  signal  a  is  recorded.  Data  in  the  first  category  may  be  divided  into  two 
types:  (i)  data  in  which  a  single  inspection  is  made  for  each  flaw,  and  (ii)  data  in  which  multiple 
inspections  are  made  for  each  flaw.  Analysis  techniques  for  estimating  POD  functions  are 
described  in  this  section.  Two  techniques  for  analyzing  NDI  results  recorded  in  the  pass/fail 
form  are  presented  first  followed  by  estimation  of  POD  functions  from  a  versus  a  data. 

2.3.1  Analysis  of  Pass/Fail  Data 

Traditionally,  NDI  reliability  data  has  been  collected  as  a  crack  length,  a„  along  with  an 
indication  of  whether  or  not  the  crack  was  found  during  a  particular  inspection.  Crack  lengths  a*, 
are  determined  through  independent  means  such  as  replicates  or  tear-down  inspections.  During 
the  inspections,  the  inspectors  record  whether  each  site  or  flaw  passed  or  failed  the  inspection. 
Because  most  of  the  NDI  reliability  data  currently  available  is  in  this  pass/fail  format,  the 
analysis  of  pass/fail  data  is  discussed  first. 

An  NDI  reliability  demonstration  experiment  can  be  conducted  in  two  ways:  one  inspection 
per  flaw  or  multiple  inspections  per  flaw.  For  data  collected  with  a  single  inspection  per  flaw, 
all  the  observations  are  independent  and  the  analysis  is  reasonably  simple.  Multiple  inspections 
conducted  on  the  same  flaw  will  be  correlated  so  that  there  are  dependencies  between  observations 
when  more  than  one  inspection  is  made  for  each  flaw.  These  two  types  of  experimental  data 
will  be  analyzed  differently. 

Two  techniques  can  be  used  to  analyze  pass/fail  data,  depending  on  the  type  of  data.  A 
regression  analysis  can  be  used  to  estimate  the  parameters  of  the  POD  model  when  there  are 
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multiple  inspections  for  each  flaw.  For  data  with  a  single  inspection  per  flaw,  the  maximum 
likelihood  method  provides  good  estimates  of  the  POD  model  parameters. 

The  analyses  described  on  the  following  are  based  on  the  log  odds  function  given  in  Eq. 
2.1a.  A  direct  analysis  of  the  model  when  expressed  in  the  form  given  by  Eq.  2.1a  is  very 
complicated.  The  analysis  can  be  simplified  by  using  the  reparameterized  model. 


POD  (a)  = 


exp  (a  +  pin  (a)) 

1  +exp(a  +  /3/n(a)) 


The  relationship  between  n  and  a  of  Eq.  2.1a  and  a  and  0  of  Eq.  2.4  is: 


n  =  -a/ 0 


=  ir  / 


For  both  the  regression  analysis  and  the  maximum  likelihood  method,  estimates  of  /i  and  a  can 
be  calculated  by  substituting  the  appropriate  estimates  of  a  and  /?  into  the  right-hand  sides  of 
Eqs.  2.5  and  2.6. 

23.1.1  Regression  Analysis 

Regression  analysis  can  be  used  for  NDI  reliability  data  in  which  (i)  multiple  inspections 
are  performed  for  each  flaw,  and  (ii)  a  single  inspection  is  performed  for  each  flaw  but  the  data 
can  be  grouped  conveniently  into  crack  length  intervals.  The  log  odds  transformation  converts 


Eq.  2.4  to 


,  POD  (a) 

'nli-p6b(a)j=a  +  ^'"‘' 


Y  =  px  +  a 

in  which  Y  (a)  and  X  are  transformed  varirbles 


»>-(r3Sfe) 
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and 


X  —  In  a 


(2.10) 


Thus,  the  linear  regression  method  can  be  used  to  estimate  a  and  0. 

Before  performing  linear  regression  on  NDI  reliability  data,  the  data  must  be  reduced  to  a 
set  of  n  pairs,  (a,,  p,)  where  ai  is  the  crack  length  for  the  i-th  pair  and  pt  is  the  proportion  of 
times  the  flaw  (or  flaws)  were  detected.  If  the  data  contain  multiple  inspections  of  each  flaw, 
ai  will  be  the  length  of  a  single  flaw  and  p,  will  be  the  proportion  of  time  that  the  flaw  was 
detected.  If  flaws  are  grouped  into  crack  length  intervals,  a,  will  be  the  midpoint  of  the  i-th 
interval  and  p,  will  be  the  proportion  of  flaws  in  the  i-th  interval  that  were  detected. 

Given  the  n  pairs  of  (a,,  p,)  data  points  to  be  fit  by  the  regression  analysis,  the  transformations 
of  Eqs.  2.9  and  2.10  are  performed,  resulting  in  a  set  of  n  (X,,  Yj)  pairs,  i.e.,  X,  = 
lna{  and  YJ  =  In  [p,  /  (1  —  pi)\ . 


Variables  X  and  Y  are  then  used  in  a  linear  regression  analysis,  Eq.  2.8,  resulting  in 
estimates  of  a  and  /?  as 


ZXiYi-nXY 

»= l 


(2.11) 


where 


d  =  Y-0X 


n  '  n  f— ' 

i=i  i=l 


(2.12) 

(2.13) 


The  estimated  mean  of  Y  as  a  function  of  a  follows  from  Eqs.  2.8  and  2.10  as 


Y  (a)  =  a  +  $  Ina 


(2.14) 


The  formula  for  a  lower  confidence  bound  on  the  mean  of  Y  (a)  is  given  by 
Yl{a)  =  a  +  0lna 


*(»-2),7(S) 


(/na  —  X)2 


SSX 


(2.15) 
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where  7  =  confidence  coefficient  (level),  /(n_2)i7  =  7-th  percentile  of  t  distribution  with  n  -  2 
degrees  of  freedom. 


(2.16) 


(2.17) 


The  inverse  transformation  of  Eq.  2.9  gives  the  estimate  of  POD  (a)  and  its  lower  confidence 
bound,  denoted  by  PODi  (a),  as  follows 


POD  (a)  = 


exp  (F  (a)) 

1  +  exp  (7  (a)) 


(2.18) 


PODL(a)= 

1+ exp  (Ft  (a)) 

in  which  Y  (a)  and  Yl(cl)  are  given  by  Eqs.  2.14  and  2.15,  respectively. 


(2.19) 


When  the  observed  proportion  of  detected  cracks  at  a  crack  length  is  zero  or  one,  i.e.. 
Pi  =  0  or  pi  =  1,  the  transformation  for  Y,  =  In  [p,/  (1  —  p,)]  is  undefined.  There  are  several 
alternatives  to  circumvent  this  problem.  One  possibility  is  to  use  the  mean  estimate,  p,  for  the 
proportion,  p,,  of  detected  cracks  at  a  crack  length. 

It  ..  .  n  i 

jr+T*  2 

5  if  ;  =  n/2  |  (2.20) 

t  +  1  .*  .  n 

rTfT  ^  <  2  i 

Another  possibility  is  to  use  the  median  estimate  p  : 

(t  -  0.3  ..  .  n 

^Tcu,f,>2 

i  if  i  =  n/2  (2.21) 

«  +  0.7  .  n  I 

l^rTo4lf,<  2 1 


(2.21) 
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Consequently,  it  follows  from  Eqs.  2.5  and  2.6  that  estimates  of  the  location  and  scale 
parameters,  //  and  a,  for  the  log  odds  model  given  by  Eq.  2.1a  are  given  by 


-&r» 

(2.22) 

1  (M) 

(2.23) 

Inspection  data  for  magnetic  particle  examinations  applied  under  water  for  detection  of 
fatigue  cracks  and  artificial  defects  was  obtained  in  graphical  form  as  shown  in  Fig.  2.4.  In 
this  figure,  each  circle  represents  a  data  point  with  the  area  of  the  circle  proportional  to  the 
number  of  observations  for  each  crack.  Unfortunately,  the  raw  data  set  was  not  available.  For 
illustrative  purposes,  the  circles  in  Fig.  2.4  were  read  graphically  without  accounting  for  the 
size  of  the  circles.  The  results  were  plotted  in  Fig.  2.5  as  circles.  Using  the  regression  analysis 
presented  above,  the  POD  (a)  function,  i.e.,  the  mean  curve,  is  plotted  as  Curve  1  and  the  lower 
95%  confidence  bound  is  plotted  as  Curve  2.  The  lower  confidence  bound.  Curve  2,  lies  above 
many  of  the  individual  data  points.  Note  that  the  confidence  bound  is  a  bound  on  the  mean  or 
POD  (a)  curve  of  the  detection  probabilities,  not  the  population  of  detection  probabilities. 

23.1.2  Maximum  Likelihood  Estimates 

In  using  the  method  of  regression  analysis,  grouping  of  NDI  observed  data  in  each  crack 
length  interval  is  required,  when  the  experiment  involves  a  single  inspection  per  flaw.  In  this 
fashion,  data  of  ( )  pairs  can  be  obtained.  Frequently,  however,  the  crack  length  for  NDI 
reliability  data  with  a  single  inspection  per  flaw  cannot  be  grouped  conveniently.  In  this  regard, 
the  method  of  maximum  likelihood  can  be  used  to  estimate  the  parameters  of  the  POD  (a) 
model  given  by  Eq.  2.1a.  With  such  an  approach,  the  parameters  are  estimated  which  maximize 
the  probability  of  obtaining  the  observed  data.  Unlike  the  regression  analysis,  the  maximum 
likelihood  estimates  do  not  require  grouping  of  data  with  a  single  inspection  per  flaw.  Instead, 
they  are  based  directly  on  the  observed  outcome  of  0  for  a  non-detection  and  1  for  a  detection. 
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To  find  the  maximum  likelihood  estimates  of  Eq.  2.4  from  a  sample  of  single  inspections 
of  n  cracks,  the  following  procedure  adopted  from  Cox  [2.15]  can  be  used.  The  maximum 
likelihood  estimates  d  and  0  of  a  and  0  satisfy  the  simultaneous  equations: 


A  A  exp(d  +  0/n(a,-)) 
i=]  ,=i  1  +  exp  (a  +  0ln  (a;)) 

A  A,n  (fl*)  exP  (“  +  &n  (“»')) 

gZ,/n(a,)  £  1  +  exp  (a  +  0ln  (a,)) 


(2.24) 


(2.25) 


in  which  the  observed  data  set  is  denoted  by  (a,,  Z,)  where  Zt  =  1  if  the  flaw  is  detected  and 
0  if  it  is  not. 


The  variances  and  covariance  of  the  estimates  d  and  0  are 

Var(d)  =  ’f-2 

(1  +exp(o  +  0ln(ai) 


(2.26) 


^(fn(Q,))2exp(a  +  ^n(qt)) 
If  (1  +exp(a  +  0/n(a;)))2 


Cov  fdj)  =  V  in  (a;)  exp  (a  +  0ln(ai)) 
'  '  (1  +  exp(a  +  0ln(ai)))2 


(2.27) 


(2.28) 


V  ’  •  ^(1+exp  {a+0ln(ai)))2 

Estimates  of  these  variances  and  covariance  are  calculated  by  substituting  the  estimates  a  and  0 
in  Eqs.  2.24  and  2.25  for  a  and  0  in  Eqs.  2.26-2.28. 

The  maximum  likelihood  estimate  of  the  POD  function  is  calculated  by  substituting  a  and  0 
for  o  and  0  in  Eq.  2.4.  The  change  of  variables  must  be  made  using  the  same  transformation 
that  was  used  in  the  regression  analysis  to  obtain 


'n[r^^)]=r(a)=“+^w 


(2.29) 


For  very  large  sample  sizes,  estimates  of  the  variances  and  covariance  of  a  and  0  can  be  used 
to  calculate  a  lower  confidence  bound  on  Y  (a)  as  given  by 


Yl(a)  =  a  +  0ln(a)  -  zyJSl  +  2ln(a)  +  5^  +  (fn(a))2Sj 


(2.30) 
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where  7  is  the  confidence  level,  z 7  satisfies  $  (z7)  =  7  with  $  ()  being  the  stan¬ 
dard  normal  distribution  function,  S\  is  the  estimate  of  Var  (a)  given  by  Eq.  2.26, 

is  the  estimate  of  Cov  ^d,  0j  given  by  Eq.  2.28,  and  S|  is  the  estimate  of  Var  given 
by  Eq.  2.27. 

Since  the  log  odds  transformation  is  monotonic,  the  inverse  transformation  of  the  confidence 
bound  on  Y  (a)  will  be  the  confidence  bound  on  POD  (a).  Specifically, 


POD  (a)  = 

1  +  exp  (K  (a)) 


(2.31) 


p°DL  (a)  =  rfP(^(°))  (2.32) 

1  +  exp  (Yl  (a)) 

in  which  Y  (a)  and  YL  (a)  are  given  by  Eqs.  2.29  and  2.30,  respectively. 

Generally,  maximum  likelihood  estimates  are  better  than  regression  estimates  that  require 
the  grouping  of  data  for  single  inspection  per  flaw;  however,  if  the  number  of  flaws  is  very  large 
(greater  than  100)  and  the  groupings  do  not  result  in  many  0’s  and  l’s  for  p,’s,  the  results  of 
both  analyses  should  be  similar. 

The  equations  for  solving  the  maximum  likelihood  estimates  d  and  0  given  by  Eqs.  2.24 
and  2.25  are  nonlinear.  In  solving  simultaneous  nonlinear  equations,  suitable  initial  estimates  of 
a  and  0  are  needed.  It  is  possible  that  d  and  0  have  more  than  one  solution  and  the  iteration 
procedures  will  converge  to  the  solution  closest  to  the  initial  estimates.  As  a  result,  the  initial 
estimates  of  d  and  0  are  important.  In  this  connection,  the  moment  method  has  been  suggested 
in  Ref.  2.7  to  determine  the  initial  estimates. 

An  example  for  the  application  of  the  method  of  maximum  likelihood  is  given  in  Fig.  2.6. 
In  this  figure,  the  circles  represent  a  set  of  hypothetical  inspection  data  for  a  single  inspection  per 
flaw.  Hence,  the  data  is  binary,  i.e.,  one  for  detection  and  zero  for  nondetection.  The  POD  (a) 
curve  is  shown  by  Curve  1  and  the  95%  lower  confidence  bound  is  indicated  by  Curve  2. 
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2.3.2  Analysis  of  a  Versus  a  Data 

The  POD  is  the  measure  of  inspection  uncertainty,  not  the  cause.  Causes  of  uncertainty 
should  be  defined  in  terms  of  the  inspection  process.  Typical  NDI  systems  apply  a  stimulus  to 
a  suspect  area  and  record  the  signal  that  returns  from  the  specimen.  A  positive  flaw  indication 
occurs  if  the  signal  is  higher  than  a  threshold  value.  The  variability  of  the  response  signals 
comes  from  the  following  reasons: 

1.  Material  variability  results  in  unpredictable  changes  in  the  stimulus  before  it  reaches  the  flaw 
and  in  the  signal  before  it  returns  to  the  NDI  system; 

2.  Variability  in  flaw  geometry  and  orientation  produces  variability  in  the  signal; 

3.  Calibration  changes  in  the  instruments  from  inspection  to  inspection  reduces  the  predictability 
of  the  signal. 

Since  a  flaw  is  detected  if  the  response  signal  a  is  larger  than  the  threshold,  the  POD  (a) 
is  the  probability  that  the  response  signal  a  is  greater  than  the  threshold  a^.  Furthermore,  the 
variability  of  the  response  signal  a  depends  on  the  inspection  process.  A  typical  plot  for  the 
response  signal  a  versus  crack  length  a  is  shown  in  Fig.  2.7.  It  is  observed  from  Fig.  2.7  that 
the  a  values  from  a  single  flaw  are  typically  grouped  around  a  point  that  is  shifted  from  the 
mean  curve.  This  pattern  of  grouping  indicates  that  there  are  two  sources  of  variation  in  the 
response  signal.  One  source  is  the  variability  in  the  mean  a  from  flaw  to  flaw,  and  the  other  is 
the  variability  in  a  from  inspection  to  inspection  of  the  same  flaw. 

The  causes  of  uncertainty  can  be  grouped  into  two  sources  of  variation:  (1)  the  material 
properties,  flaw  location,  geometry  and  orientation,  and  the  pattern  of  residual  stresses  which  are 
strictly  associated  with  individual  flaws  and  which  do  not  change  from  inspection  to  inspection; 
(2)  factors  that  change  from  inspection  to  inspection  including  human  factors,  such  as  transducer 
variability  and  calibration.  Because  of  these  two  distinct  sources  of  variation,  the  response  signal 
has  a  compound  distribution.  First  a  flaw  is  picked  at  random  along  with  its  individual  mean 
a.  Then  the  human  factors  and  equipment  factors  come  into  play  resulting  in  random  deviation 
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from  the  flaw  mean  for  an  individual  inspection.  These  are  two  distinct  random  processes  with 
distinct  random  variables.  As  a  result,  the  response  signal  a  can  be  expressed  as 


a  =  f  (a) +  c  +  e 


(2.33) 


where  /  (a)  represents  the  overall  trend  in  a  as  a  function  of  a,  c  represents  the  flaw  to  flaw 
variation,  and  e  represents  the  variation  from  inspection  to  inspection  of  the  same  flaw.  The 
function  /  (a)  is  fixed  while  variables  c  and  e  are  random  with  means  of  0. 

NDI  uncertainty  is  attributed  to  random  variation  in  the  response  signal  or  a  value  for  an 
NDI  system.  The  POD  can  be  expressed  as  the  probability  that  a  is  bigger  than  the  detection 
threshold  a*/,.  The  analysis  of  a  versus  a  data  and  estimation  of  the  POD  function  will  be 
described  in  the  following. 

Equation  2.33  provides  the  basic  model  for  the  analysis  of  a  versus  a  data.  The  flaw-related 
and  flaw-independent  terms  c  and  e  are  random  variables  with  means  equal  to  0  and  variances 
equal  to  si  and  sj,  respectively.  The  mean  and  variance  of  a  for  a  single  inspection  of  a  flaw 
of  size  a  picked  at  random  are: 


E(a/a)  =  /  (a) 
VK{a\a)  =  s2c  +  sl 


(2.34) 


Methods  for  analyzing  the  inspection  data  a  versus  a,  are  based  on  a  model  of  Eq.  2.33 
available  in  the  literature.  However,  the  choice  of  a  suitable  method  depends  on  the  functional 
form  of  /  (a),  Eq.  2.33  and  the  number  of  inspections  per  flaw.  In  general,  it  would  be  more 
convenient  to  choose  a  linear  function  for  /  (a)  through  appropriate  transformations  of  a  and  a. 
Suppose  a  and  a  are  transformed  into  Y  and  X  through: 

Y  =  Ina ,  X  =  Ina  (2.35) 

and  consider  a  linear  function  for  /  (a),  i.e.,  /  (a)  =  a  +  flX,  so  that  Eq.  2.33  becomes: 

Y  =  c*  +  0X  +  c  +  e  (2.36) 
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where  a  and  0  are  parameters  to  be  determined  from  NDI  reliability  data  of  a  versus  a.  Note 
that  a  and  0  in  Eq.  2.36  are  new  parameters  different  from  those  discussed  previously,  such 
as  those  appearing  in  Eqs.  2.4-2.6. 


If  the  error  variables  c  and  e  are  assumed  to  have  normal  distributions  with  zero  means, 
it  follows  from  Eq.  2.36  that  Y  is  a  normal  random  variable  with  mean  value  of  or  + 
0X  and  standard  deviation  S  =  (s*  +  Then,  the  POD  function  is  given  by: 

POD  (a)  =  P  (a  >  ath)  =  P  {In  (a)  >  ln(ath)}  =  P(Y>  Ytk) 


=  1 


Ytk  -  (a  +  fiX) 


(2.37) 


in  which: 


S  =  yjs\  +  a*  and  Ytk  =  ln&tk  (2.38) 

and  9  (x)  is  the  standard  normal  distribution  function.  Using  the  symmetry  properties  of  9  (x), 
Eq.  2.37  becomes: 

(x-  faf*) 

roD(°)  =  «(  S/fi 

Equation  2.39  is  a  form  of  the  lognormal  distribution  function  with  mean  jl  and  standard  deviation 
a  of  log  crack  length  given  by: 


) 


(2.39) 


H  -  (In  (ath)-  a)  IP 


(2.40) 


a  =  S/P  (2.41) 

In  the  previous  analysis,  the  log  logistic  function  of  Eq.  2.1a  was  used  to  model  the  POD  function; 
however,  the  log  logistic  function  is  a  close  approximation  to  the  lognormal  distribution.  The 
use  of  the  lognormal  distribution  above  (Eq.  2.39)  instead  of  the  log  logistic  distribution  will 
therefore  result  in  very  similar  estimates  of  the  POD  function. 

Since  Eq.  2.36  is  linear,  the  method  of  linear  regression  described  previously  can  be  used 
conveniently  to  estimate  the  parameters  a,  0  and  5  appearing  in  the  POD  (a)  function  of  Eq. 


2.39  as  follows.  The  NDI  reliability  data  is  expressed  in  n  (a,,  a,)  pairs.  These  n  pairs  of 
data  are  transformed  into  a  set  of  n  (Yi,Xi)  pairs  through  the  transformation  of  Eq.  2.35, 
i.e.,  Y,  =  Inhi  and  X{  =  /na,.  By  use  of  linear  regression  analysis,  formulas  for  estimates 
a,  P  and  S  for  parameters  a,  p  and  5  are  identical  to  those  given  by  Eqs.  2.11,  2.12  and 
2.16,  respectively.  After  estimating  a,  p  and  S  from  these  equations,  parameters  ji  and  «r 
can  be  computed  from  Eqs.  2.40  and  2.41  in  which  a,  P  and  5  are  replaced  by  a,  p  and  S, 
respectively.  Hence,  the  POD(a)  function  (Eq.  139)  can  be  expressed  as: 

*»(„)  =  (2.42) 


It  is  noted  from  Eqs.  2.39-2.42  that  the  POD(a)  function  depends  on  the  specifying 
detection  threshold  a*&.  The  effects  of  the  detection  threshold  at *  on  POD(a)  are  described 
in  the  following.  First,  the  median  detection  crack  length  increases  with  the  detection  threshold. 
Second,  the  slope  of  the  POD(a)  function  decreases  as  the  detection  threshold  increases. 

A  method  described  by  Cheng  and  lies  [2.16]  for  calculating  confidence  bounds  on  die 
lognormal  cumulative  distribution  can  be  adopted  to  calculate  confidence  bounds  on  the  POD 
function  given  by  Eq.  2.42.  The  formulas  given  by  Cheng  and  lies  [2.16]  for  the  7  percent 
lower  confidence  bound  can  be  used  as  follows: 


POD(a)  =  $(*i) 


<2. 43) 


in  which 


I  A 

=w 


+1£zZ)!+1 

+  SSX  * 


) 


(2.44) 


where 


__  1  » 

X  =  Ina  and  X  =  —  Inai  (2.45) 

isl 

In  Eq.  2.44,  n  is  the  sample  size,  A  is  the  7  -  th  percentile  of  a  Chi-Square  distribution  with 


two-degrees-of-freedom,  SSX  is  given  by  Eq.  2.17  and: 

X-fi 


z  = 


(2.46) 
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Theoretically  speaking  the  a.  versus  a  data  sets  are  superior  to  the  pass/fail  data  sets,  since 
a  versus  a  data  contains  more  statistical  information  regarding  the  uncertainties  of  an  NDI 
system.  With  the  current  state  of  NDI  technology  for  marine  structures,  extensive  research 
effort  is  needed  to  obtain  a  versus  a  data. 

2.4  Conclusions 

A  review  of  POD  curves  used  in  the  aerospace  industry  was  presented  in  this  chapter.  This 
review  was  done  because  it  is  believed  that  the  shape  of  POD  curves  used  in  the  aerospace 
industry  can  provide  useful  guidelines  for  assessing  the  reliability  of  flaw  detection  and  for 
establishing  POD  curves  for  marine  structures.  Emphasis  was  given  to  the  log  odds  model 
which  was  investigated  extensively.  Another  very  good  model  is  the  Weibull  model  which  has 
been  well  studied  in  dealing  with  other  applications  and  found  to  be  very  reasonable  for  die 
POD  function.  Indeed,  it  is  a  special  case  of  the  Weibull  model  that  will  be  used  as  a  POD 
curve  in  the  third  chapter  that  follows. 
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PROBABILITY  OF  DETECTION 


0 


Fig.  2.3  Schematic  of  Probability  Density  Function  of  Crack  De¬ 
tection  Probabilities  at  a  Crack  Length. 


Ref:  Shell  ♦  Kn 


Fig.  2.4  Reliability  of  Magnetic  Particle  Examination  Applied  Un¬ 
der  Water  for  Detection  of  Fatigue  Cracks  and  Artificial 
Defects.  Observation  Points  are  Given  with  an  Area  Pro¬ 
portional  to  the  Number  of  Individual  Observations  Con¬ 
tained  in  Each  Point. 
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Curve  1 
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Fig.  2.5  Regression  Estimate  of  POD  Function  for  Pass/Fail  Data 
with  Multiple  Inspections. 
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Fig.  2.6  Maximum  Likelihood  Estimate  of  POD  Function  for 
Pass/Fail  Data  with  a  Single  Inspection. 
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Fig.  2.7  Observed  Peak  Eddy  Current  Response  Voltages  From 
Two  Measurements  on  Eack  Crack. 
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III.  STRUCTURAL  RELIABILITY  UNDER 
BAYESIAN  INSPECTION 


3.1  Introduction 

After  the  first  two  introductory  chapters,  the  third  chapter  constitutes  the  main  part  of 
this  work  whose  basic  objective  is  to  develop  a  non-periodic  inspection  procedure  for  marine 
structures  so  as  to  maintain  their  reliability  at  a  prespecified  design  level  throughout  their  lifetime. 
In  this  work,  the  main  cause  of  a  structure’s  reliability  decline  between  inspections  is  considered 
to  be  the  formation  and  propagation  of  fatigue  cracks.  Fatigue  is  actually  one  of  the  most 
important  problems  of  offshore  structures  located  in  harsh  environments  such  as  the  North  Sea. 
The  problem  of  fatigue  is  further  complicated  by  the  random  nature  of  the  loading  to  which 
these  structures  are  subjected. 

Fatigue  damage  is  considered  to  be  initiated  in  a  structure  when  the  smallest  size  measurable 
crack  develops,  whether  or  not  it  is  detected.  The  fatigue  process  in  a  structural  member  consists 
of  crack  initiation,  followed  by  crack  propagation  and  the  resulting  member  strength  degradation 
(see  Fig.  3.1).  Periodic  inspections  of  fatigue-sensitive  structures  have  been  common  practice  in 
order  to  maintain  the  reliability  of  the  structures  at  the  desired  prespecified  level;  if  a  fatigue  crack 
is  detected  by  inspection,  the  cracked  component  is  repaired  or  replaced  with  a  new  member  so 
that  both  the  residual  strength  and  fatigue  characteristics  of  this  component  are  renewed,  thus 
increasing  the  overall  reliability  of  the  structure.  Hence,  reliability  analysis  of  fatigue-sensitive 
structures  under  random  loading  and  periodic  inspections  is  of  practical  importance  and  forms 
the  primary  concern  of  this  study. 

In  previous  studies  performed,  for  example,  by  Yang  and  Trapp  [3.3]  and  Paliou  and 
Shinozuka  [3.2],  the  effect  of  periodic  inspections  on  fatigue-sensitive  structures  was  examined 
with  probabilistic  methods.  These  studies  paved  the  way  for  the  development  of  reliability-based 
design  criteria  for  ocean  structures  involving  the  notion  of  periodic  or  nonperiodic  inspections. 
However,  further  study  is  needed  to  deal  with  (a)  the  lack  of  inspection  data  with  which  statistical 
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analyses  can  be  performed,  (b)  uncertainty  in  the  fatigue  properties  of  structural  components,  (c) 
difficulties  in  estimating  the  stresses  acting  on  members  due  to  structural  complexity,  (d)  errors 
during  construction,  and  (e)  uncertainty  in  the  detection  capabilities  of  the  various  inspection 
techniques  to  be  used.  In  fact,  Itagaki  and  Yamamoto  [3.1]  analyzed  the  effect  of  non-periodic 
inspections  using  Bayesian  analysis  and  determined  the  inspection  interval  so  as  to  maintain  the 
reliability  of  the  whole  structure  at  some  prespecified  design  level.  The  present  study  represents 
a  further  extension  of  this  work  by  Itagaki  and  Yamamoto  in  order  to  include  such  new  features 
as  the  detection  of  a  crack  of  a  certain  length.  Monte  Carlo  simulations  are  carried  out  to  verify 
the  validity  of  the  methodology. 

Section  3.2  describes  the  analytical  models  used  and  the  various  assumptions  employed  in 
this  analyse,  including  expressions  for  fatigue  crack  initiation  time,  fatigue  crack  propagation  and 
failure  rate.  In  Section  3.3,  analytical  interpretation  of  inspection  results  involving  definition  of 
all  possible  outcomes  of  such  inspections  is  given.  According  to  these  definitions,  the  probability 
of  occurrence  of  each  of  these  outcomes  and  related  events  is  presented  in  detail.  Sections  3.4 
and  3.5  are  concerned  with  the  determination  of  anpropriate  inspection  intervals  so  that  the 
structural  reliability  is  kept  at  the  desired  level.  For  this  purpose,  a  Bayesian  approach  is  applied 
to  treat  the  various  uncertainties  involved.  The  types  of  uncertainty  considered  in  this  study 
are  those  appearing  in  the  expressions  for  (a)  fatigue  crack  initiation  time,  (b)  fatigue  crack 
propagation  rate,  and  (c)  probability  of  crack  detection. 

Assuming  prior  density  functions  for  the  unknown  parameters,  the  inspection  results  are  used 
in  accordance  with  Bayes  Theorem  to  upgrade  the  prior  density  functions.  A  general  formulation 
is  given  for  the  case  where  a  detailed  record  of  the  entire  inspection  history,  including  repair 
or  replacement  records  for  each  and  every  member,  is  available.  Numerical  simulations  were 
carried  out  (see  Section  3.6)  to  verify  the  validity  of  this  Bayesian  approach 

3.2  Basic  Assumptions 

The  following  assumptions  are  made  for  the  purposes  of  this  study: 
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1.  In  each  structural  member  there  is  only  one  fatigue-critical  location  (hotspot)  where  a 
crack  can  initiate. 

2.  All  structural  members  are  inspected  immediately  after  initiation  of  service  and  at  the 
time  of  each  scheduled  inspection.  If  a  member  is  found  not  to  be  intact,  the  following 
action  is  taken: 

•  If  a  crack  is  detected  in  a  member,  that  member  is  repaired  and  regains  its  initial 
strength  characteristics. 

•  If  the  member  is  found  to  have  failed,  it  is  replaced  by  a  new  one. 

3.  The  entire  inspection  history  of  each  member  is  considered  to  be  known  at  the  time  of 
the  current  inspection. 

4.  For  fatigue  crack  initiation,  the  time  to  crack  initiation  (TTCI),  denoted  by  t,  is  assumed 
to  be  a  random  variable  with  a  density  function  following  the  Weibull  distribution: 

i>o  <3i> 

The  uncertainty  in  the  TTCI  is  introduced  by  the  scale  parameter  /?.  Hence,  Eq.  3.1 
indicates  a  Weibull  density  conditional  to  a  given  value  of  /?.  The  shape  parameter  a 
is  assumed  to  be  deterministic  for  the  sake  of  simplicity.  The  distribution  function  of 
the  TTCI  is  expressed  by: 

Fc(t  1 0)  =  1  -exp  f(^j  1  <>0  (3.2) 

5.  For  fatigue  crack  propagation,  fracture  mechanics  theory  is  used  to  determine  the  length 
of  a  propagating  crack  under  random  stress.  It  is  assumed  that  the  crack  grows  in 
accordance  with  the  following  law: 

Yt  =  c\  (Atf)2  =  Ci  (A<r  •  v/^)2  =  ca  (3.3) 
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with 


c  =  irCiAa2  (3.4) 

where  a  is  the  crack  length,  C\  is  a  material  constant,  A K  is  the  stress-intensity-factor 
fluctuation  and  Aa  is  the  nominal  stress  fluctuation.  Integrating  Eq.  3.3  from  the  initial 
crack  length  ao  at  the  TTCI  =  tc,  up  to  the  current  crack  length  a(t  -  tc)  at  time  t,  the 
following  result  is  obtained: 

a{t  -tc\c)  -  a0exp  [c(f  -  fc)]  (3.5) 

The  uncertainty  in  fatigue  crack  propagation  is  introduced  by  parameter  c,  and  therefore 
the  crack  length  indicated  by  Eq.  3.5  is  conditional  to  a  given  value  of  c. 

6.  The  probability  of  detecting  a  fatigue  crack  of  length  a  at  the  time  of  inspection  is  given 
by: 

D  (a  |  d)  =  1  —  exp  [— d(a  —  ao)]  (3.6) 

The  uncertainty  in  the  probability  of  crack  detection  is  introduced  by  parameter  d.  Thus, 
the  probability  of  detection  shown  by  Eq.  3.6  is  conditional  to  a  given  value  of  d.  The 
minimum  detectable  crack  length  is  denoted  by  ao. 

7.  If  a  crack  is  detected  in  a  member  at  the  time  of  inspection,  the  crack  length  is  assumed 
to  be  accurately  measured. 

8.  Failure  of  a  member  occurs  when  random  stress  exceeds  the  strength  of  the  member  for 
the  first  time.  A  member  can  fail  either  before  or  after  crack  initiation.  With  the  former, 
the  failure  rate  is  a  constant  depending  only  on  the  characteristics  of  the  stationary 
random  process  representing  the  stress  fluctuation  at  the  hotspot  However,  with  the 
latter,  the  failure  rate  also  depends  on  the  crack  size,  cm  which  the  member’s  residual 
strength  depends.  Therefore,  the  problem  of  evaluating  the  failure  rate  is  essentially  that 
of  estimating  the  first-passage  failure  probability  with  a  constant  or  variable  two-sided 
threshold  depending  on  whether  the  member  fails  before  or  after  crack  initiation.  The 
failure  rates  before  and  after  crack  initiation  at  tc  take  the  following  forms: 
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Before  crack  initiation: 


k  ( t )  =  ho  =  exp  (r) 


(3.7) 


After  crack  initiation: 


h  ( t )  =  exp  [g  (t  -  tc)  +  r]  (3.8) 

For  the  sake  of  simplicity,  parameters  r  and  g  are  assumed  to  be  deterministic. 

The  reliability  of  a  member  before  crack  initiation  during  the  service  period  from  I)  to 
t  is  denoted  by  U  (t  -  T/)  and  given  by: 


U(t  —  Tj)  =  exp  < 


-  J  =exp|-  J  hodrj  =exp|-  yexp(r)dr| 

(3.9a) 


or 


U  (t  -  Ti)  =  exp  {-  (t  -  Ti)  •  exp  (r)}  for  t<tc  (3.9b) 

where  7)  is  the  time  of  service  initiation  for  the  member  under  consideration;  this  implies 
that  the  member  was  repaired  or  replaced  at  the  time  of  the  l — th  inspection. 

The  reliability  of  a  member  after  crack  initiation  during  the  service  period  from  te  to  t 
is  denoted  by  V  (t  -  tc)  and  given  by: 


V(t 


-tc)  =  cxpl^-  J  h(r)dT  =  exp  j  -  ^ 


exp  [q(r  -  tc)  +  r)dr 


}  0 


10a) 


or 


V(t-  tc)  =  exp  [exp  {q  ( t  -  tc)  +  r}  -  exp  (r}]|  for  t>tc  (3.10b) 


9.  The  probability  of  detecting  member  failure  at  the  time  of  inspection  is  equal  to  one, 
if  such  a  failure  indeed  exists. 

10.  No  stress  redistribution  is  considered  in  the  structure  after  the  occurrence  of  member 
failure  or  during  crack  propagation. 
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3.3  Possible  Events  at  Time  of  Inspection 


33.1  Definitions 

At  the  time  of  the  j -th  inspection,  Ty,  of  a  certain  member  (with  the  knowledge  that  this 
member  was  repaired  or  replaced  immediately  as  a  result  of  the  1-th  inspection  performed  at 
time  T;  (/  <  j  -  1)  or  that  this  member  initiated  service  at  time  7}  denoting  the  beginning  of 
service  for  the  structure),  one  of  the  following  three  events  may  occur 

1.  {A  :  j ,  /}  =  event  that  the  member  is  found  to  have  failed  at  the  time  of  the  j— th 
inspection  Ty  or  equivalently,  the  event  that  failure  of  the  member  occurs  during  the 
time  interval  [Ty_i,Ty]. 

This  event  consists  of  the  following  two  mutually  exclusive  events: 

Eij  =  event  that  the  member  fails  before  crack  initiation,  sometime  during  the  time 
interval  between  the  two  consecutive  inspections  at  Ty_ i  and  Ty 
Eiyj  —  event  that  the  member  fails  after  crack  initiation,  sometime  during  the  time 
interval  between  the  two  consecutive  inspections  at  T}-\  and  Tj. 

It  is  assumed  that  although  member  failure  is  always  detectable  at  die  time  of  inspection 
if  such  a  failure  indeed  exists,  it  is  impossible  to  determine  to  which  of  these  two 
events,  E\,j  and  fyj,  the  failure  Delongs. 

2.  {B\  (ctj) :  j,  1}  —  event  that  the  member  is  found  not  to  have  failed  at  the  time  of  die 
j — th  inspection  Tj  and  a  crack  of  length  between  ctj  and  ctj  +  dctj  is  detected  in  the 
member. 

Event  {B\  (ay) :  j,l}  is  alternatively  denoted  by  Ezj. 

3.  {2? 2 :  j,  /}  =  event  that  the  member  is  found  not  to  have  failed  at  the  time  of  the  j — th 
inspection  Ty  and  no  crack  is  detected  in  the  member. 

This  event  consists  of  the  following  two  mutually  exclusive  events: 

E\  j  =  event  that  the  member  does  not  fail  in  the  time  interval  [3y_i,  Ty)  and  no  crack 
exists  in  the  member  at  the  time  of  inspection  Ty. 
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E$j  =  event  that  the  member  does  not  fail  in  the  time  interval  [2j_i,Tj]  and  a  crack 
exists  in  the  member  which  is  not  detected  at  the  time  of  inspection  Ty 

In  the  following,  the  probabilities  of  these  events  for  a  particular  member  will  be  evaluated  in 
terms  of  the  probability  density  and  distribudon  functions  fc  (t  |  0)  and  Fc  (t  |  0)  of  the  TTG, 
reliability  functions  U  ( t )  and  V  (t)  and  probability  of  crack  detection  D  (a  )  d  |  c),  with  the 
knowledge  that  this  member  was  repaired  or  replaced  as  a  result  of  the  l-th  inspection  performed 
at  Ti  (l  <  j  —  1)  or  that  this  member  initiated  service  at  time  7}  denoting  the  beginning  of  service 
for  the  structure. 

3.3.2  Evaluation  of  Probabilities  of  Various  Events 

33.2.1  Event  Eij 

Event  Eij  consists  of  two  mutually  exclusive  events,  E\  -  and  E\j,  which  are  defined  as  follows: 

1.  Efj  =  event  that  a  crack  initiates  after  Tv  and  the  member  fails  before  crack  initiation 
sometime  during  the  time  interval  [T}_i,Tj].  The  probability  Pfj  of  event  E\}  is  given 
by: 

pi,j  =  {1  -Fc  (T,-  -  T,  |  0)}  ■  {V  (Ti-i  -T,)-U  (Tj  -  T,))  (3.11) 

where  {1  -  Fc  (T}  -  Tj  |  0)}  denotes  the  probability  that  the  crack  will  initiate  af¬ 
ter  Tj  and  {17  (T,_i  -  Ti)  -  U  ( Tj  -  7j)}  denotes  the  probability  that  the  mem¬ 
ber  will  fail  sometime  during  the  time  interval  [Tj_i,7,J].  The  latter  is  the 
conditional  probability  given  that  a  crack  initiates  after  T},  and  is  found  from 
{1  -  U(Tj  -  T,)  -  [1  -U(Tj-i  -  Ti)]}  (see  Fig.  3.2a). 

2.  E\ j  =  event  that  a  crack  initiates  at  some  time  instant  t  in  the  time  interval  [T;_i,T,] 
and  the  member  fails  before  crack  initiation  sometime  during  the  interval  [Tj_i,t]. 

The  probability  Ki  of  event  E*  }  is  given  by: 

T, 

Pfj=  J  fc(t  -  Ti  \  /))  ■  [U  (T,-,  -  T,)  -  U  (t  -  r,))  it  (3.12) 

T,-, 
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where  fc  (<  —  T>  |  0)  dt  denotes  the  probability  that  a  crack  will  initiate  during  the  time 
interval  [<,<  +  dt]  and  {27  {Tj-\  -Tfi-U  (t-  7})}  denotes  the  probability  that  the 
member  will  fail  during  the  time  interval  (see  Fig.  3.2b). 


Finally,  the  probability  P\<}  of  event  E\j  is  given  by: 

Pu  =  f?j  +rtj  = 

=  {1  -Fc (T,  - T,  |  /)))  ■  {V (Tj-i  -T,)-U (Tj  - T,)} 

t,  (313) 

+  J  fc{t-T,\0)-{U{Tj-l-T,)-U{t-Tl)}dt 

tU 


33.2 2  Event  E2J 

1.  Event  E2j  consists  of  two  mutually  exclusive  events,  and  Ejj,  defined  in  the 
following: 

E\}  =  event  that  a  crack  initiates  at  some  time  instant  t  in  the  time  interval 
[Ti,Ti+]]  (t  =  /, j  -  2)  and  the  crack  is  not  detected  during  all  the  subsequent  in¬ 
spections  (from  inspection  at  time  Tj+i  up  to  inspection  at  time  Tj- 1  inclusive)  and  die 
member  fails  sometime  during  the  time  interval  [Tj-i,Tj]. 

The  probability  P of  event  E%}  is  given  by: 

i- 2  (  Titl 

X  Ti  0.14) 

•  n  {1  - - <  | c)  | <0}j  <«} 

.*=i+l  J  ) 

where  fe{t  -  T\  |  0)dt  denotes  die  probability  that  a  crack  will  initiate  during  die 
time  interval  [t,  t  +  dt],  U  (t  -  T\)  denotes  the  probability  that  the  member  will  sur¬ 
vive  during  the  time  interval  [l},*],  {V  {Tj- j  - 1)  -  V  {Tj  -  <)}  denotes  the  prob¬ 
ability  that  the  member  will  fail  sometime  during  the  time  interval  [Tj-i,Tj]  and 
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{1  -  D(a(Tk  - 1 1  c)  |  <f)}  denotes  the  probability  that  a  crack  will  not  be  detected 
at  inspection  7*  (see  Fig.  3.2c). 

2.  E\  -  =  event  that  a  crack  initiates  at  some  time  instant  tin  the  time  interval  [Tj_i ,  Tj) 
and  the  member  fails  sometime  during  the  time  interval  [<,  Tj], 

The  probability  pij  of  event  E\j  is  given  by: 

I  fc(t-T,\p)  V(t-T,)  [l-V(.Tj-t)]dt  (3.15) 

n-, 

where  fc  (t  -  T/ 1  0)  dt  denotes  the  probability  that  a  crack  will  initiate  during  the  time 
interval  [t,  t  +  dt),  U(t-  Tj)  denotes  the  probability  that  the  member  will  survive  during 
the  time  interval  [Tj,  t]  and  [1  -  V  (Tj  - 1)\  denotes  the  probability  that  the  member  will 
fail  sometime  during  the  time  interval  [t,  Tj]  (see  Fig.  3.2d). 

Finally,  the  probability  P2j  of  event  fyj  is  given  by 

P-u  =  Ptj  +  ^  = 

J  (  T‘i' 

=  £{  /  /.  (I  -  T,  |  0)  ■  (/  (<  -  T,)  •  |V  (T/_,  _  t)  _  V  (Tj  - 1)] 

«  U 

i-i 

n  u-d(<.(t,-(|<:)i<o}  a 

*•+1 

+  f  M*  -  T)  1 0)  •  17  (t  -  T,)  •  [1  -  V  (Tj  - 1))  it 

tU 

3323  Event 

Event  Eij  is  defined  as  follows: 

=  event  that  the  member  does  not  fail,  but  a  crack  of  length  between  a,  and  aj+daj 
is  detected  at  die  time  of  the  j-th  inspection,  Tj.  Since  a  crack  with  length  between 


54 


aj  and  aj  +  da}  is  found  at  the  time  of  the  j-th  inspection,  the  time  te  of  initiation  of 
this  crack  can  be  computed  from 

a  {Tj  -  <c  I  c)  =  aj  =  a0  •  exp  [c  ( Tj  -  <«)]  (3.17) 

as 

tc  =  Tj  -  i  In  (3.18) 


Then,  the  probability  jnjdaj  of  event  E$j  is  given  by: 

P3,jdaj  =  p3j  (aj)daj  =  fc  ( tc  ~  7/  |  P)  dtc  ‘U  {tc  —  3))  •  V  {Tj  —  te) 

■  [  n  {l-<  C(o(T,-(o|c)|<i)) 

U=/+ 1 

On  the  right-hand  side  of  Eq.  3.19,  tc  is  to  be  replaced  by  the  expression  shown  in  Eq.  3.18. 


D{aj\c\d) 


(3.19) 


Similarly,  on  the  right-hand  side  of  Eq.  3.19,  dte  =  \dtc/daj\ da}  =  da}/  {caj).  In  this  way,  die 

right-hand  side  of  Eq.  3.19  is  given  completely  as  a  function  of  aj.  In  Eq.  3.19, 6  is  given  by 

(  1  for  T*  >  tc  'I 

6=  \  [  (3.20) 

[  0  for  7*  <  tc  ) 

and  fc  (te  —  T|  |  0)  dtc  denotes  the  probability  that  a  crack  will  initiate  during  the  time  interval 
[<«>  tc  +  dtc].  U  (tc  -  Tj)  denotes  the  probability  that  a  member  will  survive  during  the  time 
interval  [T/,te],  V  ( Tj  -  tc)  denotes  the  probability  that  a  member  will  survive  during  the  time 
interval  [<e,  Tj],  {l-6D{a{Tk-tc\c)\ d)}  denotes  the  probability  that  a  crack  will  not  be 
detected  at  inspection  7*  and  D  {aj  |  c  |  d)  denotes  die  fnobability  that  a  crack  will  be  delected 


at  inspection  Tj  (see  Fig.  3.2e). 


3314  Event  E<j 


Event  Eij  is  defined  as  follows: 

Eij  =  event  that  member  does  not  fail  and  no  crack  exists  in  member  at  time  of 
inspection  Tj. 
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J  'I  ' 


The  probability  P4j  of  event  E4tj  is  given  by: 


P4J  =  {1  -  Fe  (Tj  -  Tt  |  ff)}  •  U  (Tj  -  T|) 


(3.21) 


where  {1  -  Fc(Tj  -  Tt  |  0)}  denotes  the  probability  that  a  crack  will  initiate  after  Tj  and 
U  (Tj  -  Ti )  denotes  the  probability  that  a  member  will  survive  the  time  interval  [Tt,Tj]. 

33.2.5  Event  Ejj 


Event  E$j  is  defined  as  follows: 

E$j  =  event  that  member  does  not  fail,  and  a  crack  exists  in  the  member  which  is  not 
detected  at  the  time  of  inspection  Tj. 

The  probability  P5iJ  of  event  E$tj  is  given  by: 
i-1  f  Tit' 

^  =  /  fc(t-Tl\P)’U(t-T,)-V(Tj-t) 

,ss/  l  i  (3.22) 


{l-D(o(2l-«|c)|d)} 


where  fc(t  -  T/ 1  0)dt  denotes  the  probability  that  a  crack  will  initiate  during  time  interval 
[M  +  dt\,  U(t~  Ti )  denotes  the  probability  that  a  member  will  survive  during  time  interval 
[T|,t],  V  (Tj  —  t)  denotes  the  probability  that  a  member  will  survive  during  time  interval  [f,  Tj), 
and  (1  -  D  (a  (T*  - 1 1  c)  |  d)}  denotes  the  probability  that  a  crack  will  not  be  detected  at 
inspection  7*  (see  Fig.  3.2f). 

Finally,  the  probabilities  of  events  {A  :  j,  l),  {Bj  (aj)  :  j,  l }  &  i  {Bj  :  j,  l }  are  given  by 

P  {B\  (aj)  :  j,l)  =  Pijdaj 
P{B2:j,l}=P4j  +  Pij 


(3.23) 


At  this  point,  is  should  be  noted  that  Pm  {A  :  j,  /},  Pm  {B\  (aj) :  j,  1}  and  Pm  {B2  :  j,  /}  are 
written  for  P  {A  :  j,  /},  P  {Bx  (a}) :  j ,  1}  and  P  {B2  : /} ,  respectively,  in  order  to  identify 
member  number  m. 

3.4  Reliability  of  Member  at  Time  Instant  V  After  i-th  Inspection 

We  distinguish  between  the  reliability  of  two  types  of  members  at  time  instant  t*  after  the 
j-th  inspection  but  before  the  ( j  + 1)  -th  inspection  (Tj  <i*  <  Tj+ 1). 

3.4.1  Members  Repaired  at  j-th  Inspection 

It  is  known  that  members  are  replaced  or  repaired  at  the  j-th  inspection  in  the 
case  of  event  {A  :  j,  1}  or  event  {Bi  (aj) :  j,l}.  Writing  R  (t* ;  Repair)  instead  of 
R  (**;  Replacement  or  Repair)  for  brevity,  the  reliability  R  (t*;  Repair)  of  a  member  of  this  type 
is  given  by  the  sum  of  the  following  two  probabilities:  (a)  the  probability  that  a  member  will 
not  fail  during  time  interval  [Tj,  <*]  and  a  crack  will  initiate  after  <*;  and  (b)  the  probability  that 
a  crack  will  initiate  during  time  interval  [Tj,  <•],  but  the  member  will  not  fail  during  the  same 
time  interval.  R(t* :  Repair)  is  calculated  as  follows: 

fl(<* :  Repair)  =  {1  -  Fe  (t*  -  Tj  |  0)}  U (f  -  Tj) 

i* 

+  j  -t)dt  0M) 

T, 

where  {1  -  Fc(t*  -  Tj  \  p)}  denotes  the  probability  that  a  crack  will  initiate  af¬ 
ter  <*,  V  {<*  -  Tj}  denotes  the  probability  that  a  member  will  survive  time  interval 
-  Tj  I  $)dt  denotes  the  probability  that  a  crack  will  initiate  during  time  inter¬ 
val  [*,<  +  «ft],  U(i-  Tj)  denotes  the  probability  that  a  member  will  survive  time  interval 
[Tj,<]  and  V(t*  - 1)  denotes  the  probability  that  a  member  will  survive  during  time  interval 
[Ml- 

3.4 2  Members  Not  Repaired  at  j-th  Inspection 

It  is  known  that  members  are  not  repaired  at  the  ;-th  inspection  in  the  case  of  event 

{^2  :j,0- 
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The  reliability  R(t* :  No  repair)  of  a  member  of  this  type  is  given  by  the  sum  of  the 
following  three  probabilities  (written  as  Z)  divided  by  (P^j  +  Psj)  which  represents  die 
probability  of  event  {B2  :  j,  l}: 

a.  Probability  that  member  will  not  fail  during  time  interval  [7),  <*]  and  crack  will  initiate 
after  t*. 

b.  Probability  that  a  crack  will  initiate  during  time  interval  [Tj,t*],  but  member  will  not 
fail  during  time  interval 

c.  Probability  that  crack  initiates  at  some  time  instant  t  during  time  interval 
[Tj,Ti+1]  (t  =  l,  ...,j  -  1)  and  this  crack  is  not  detected  during  all  subsequent  inspec¬ 
tions  (from  inspection  at  time  Tj+ \  to  inspection  at  time  Tj  inclusive)  and  member  will 
not  fail  during  time  interval 


Hence, 

R(f  ■  No  repair)  =  -f  p  0.25) 

In  Eq.  3.25, 

Z  =  {l-Fc(r-r/|^)}.B(t*-T/) 

f 

+  f  f'(t-T,\p)V(t-T,)V(t'-t)it 

T, 

+  ±{j  {1  —  D(a(T,  —  1 1  c)  |  d)}J  <ft| 

(3.26) 

where  {1  -  JPC(<*  -  Ti  \  (})}  denotes  the  probability  that  a  crack  will  initiate  after  t*,  U (t*  -  T\ ) 
denotes  the  probability  that  a  member  will  survive  time  interval  (7),  <*] ,  fe  (t  -  7)  |  0)  dt  denotes 
the  probability  that  a  crack  will  initiate  during  time  interval  (< ,  t  +  dt] ,  U  (t  -  7})  denotes  the 
probability  that  a  member  will  survive  time  interval  [T/,f] ,  V  (f*  - 1)  denotes  the  probability 
that  a  member  will  survive  time  interval  [<,*♦]  and  {1  -  D(a(Tk  - i  |  c)  |  d)}  denotes  the 
probability  that  a  crack  will  not  be  detected  at  inspection  Ti-  Note  that  R(t*; No  Repair) 
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indicates  the  probability  of  event  A  that  a  member  survives  time  interval  given  event  B 
that  it  has  not  been  replaced  or  repaired  at  die  j-th  inspection.  Therefore,  P  {.A  n  B)  =  P  {A}- 
P{B  \  A}  where  P  {A  D  b]  —  Z,  P  {A}  =  P^j  +  P$j  and  P  {B  \  A)  =  R(t*; No  Repair). 
Equation  3.25  immediately  follows  from  this. 


3.5  Bayesian  Analysis 

3.5.1  Uncertain  Parameters  and  Their  Prior  Density  Function 


In  the  present  study,  /?,  d  and  c  are  considered  as  possible  sources  of  the  uncertainty. 
Initially,  a  uniform  distribution  is  assumed  for  the  three  uncertain  parameters  having  the  following 
jointly  uniform  density  function: 


/°(M,c)  = 


where: 


(An  ax  ~  An  in)  (Anax  -  An  in)  (cmax  - 
Anin  ^  0  ^  An  ax 


=  constant 


(3.27) 


Anin  ^  d  <  Anax 


(3.28) 


^  in  ^  c  ^ 


3.5.2  Likelihood  Function  as  Result  of  j-th  Inspection 

The  likelihood  function  LFj  for  the  entire  structure  as  a  result  of  the  j-th  inspection  is 
expressed  as  follows: 

M 

= n  LF\m)  (3.29) 

ms] 

where  is  the  likelihood  function  as  a  result  of  the  j-th  inspection  for  member  m  and  M 
is  the  total  number  of  members  in  the  structure. 

For  a  typical  member  m,  assume  that  replacement  due  to  failure  or  repair  due  to  a  detected 
crack  occurred  at  the  time  of  inspections  T|nT|a,...,T|r  where  r  indicates  the  number  of  times 
the  member  has  been  repaired  or  replaced  before  the  j-th  inspection,  and: 

h<h<  ...  <lT<)  (3.30) 
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It  is  pointed  out  that  h,  /2, lr  are  all  known  at  the  time  of  the  j-th  inspection  since  the  whole 
inspection  history  of  each  member  is  considered  to  be  known.  Note  also  that: 

lr<j- 1  (3.31) 


Then,  the  likelihood  function  as  a  result  of  the  j  -th  inspection  for  member  m  is  given  by: 

r 

LFjm>  =Pm{X  :  j,  lr)  ■  n  Pm  {r  :  k,  It- 1}  (3-32) 

k=l 

It  is  noted  that  /1?  /2, lr  as  well  as  r  usually  take  certain  values  unique  to  each  member. 

In  Eq.  3.32,  X  stands  for  either  A  or  Bx  (ay)  or  B2  depending  on  the  result  of  the  ;  -th 
inspection  for  member  m.  Specifically,  if  at  the  time  of  the  j— th  inspection,  member  m  is  found 
to  have  failed,  then  X  stands  for  A;  if  member  m  is  found  to  have  a  crack  of  length  between 
ay  and  ay  +  day,  then  A'  stands  for  B\  (ay)  and  if  member  m  is  found  intact,  then  ATstands  for 
B2.  Also  in  Eq.  3.32,  Y  stands  for  either  A  or  B\  ( aik )  depending  on  the  result  of  the  /* — th 
inspection  for  member  m.  Specifically,  if  at  the  time  of  the  /* — th  inspection,  member  m  is 
found  to  have  failed,  then  Y  stands  for  A  and  if  member  m  is  found  to  have  a  crack  of  length 
between  a\k  and  aik  -f  da;fc,  then  Y  stands  for  B\  (ajk).  Finally,  for  the  case  where  member  m 
is  found  intact  at  all  inspections  prior  to  the  j— th,  the  product  appearing  in  Eq.  3.32  is  set  equal 
to  1  and  Eq.  3.32  takes  the  form: 


=Pm(X  :  j,l 0} 


(3.33) 


Note  that  Iq  denotes  the  time  of  initiation  of  service  for  the  structure. 


3.53  Posterior  Joint  Density  Function  of  Uncertain  Parameters 


The  posterior  joint  density  function  of  the  three  uncertain  parameters  immediately  after  the 
j— th  inspection  is  given  by: 


fj(P,d,c)  = 


_ LFjf0  (0,d,c) _ 

ftn  ax  4nax  <max 

Iff  (Numerator)  d0d  (d)  dc 
^min  ^nin  ‘win 


(3.34) 
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3.5.4  Reliability  of  Entire  Structure  at  Time  t* 

The  reliability  of  the  entire  structure  consisting  of  M  members  at  time  t*  after  the  j'-th 
inspection  is  denoted  by  Rm  (<*)  and  is  given  by: 

£max  dmax  nnax 

T&(t')  =  j  J  J  Rkt(t’\P^i)f‘(0,c,d)d0d(i)dc  (3.35) 


8min  ^min  cmin 


where: 


(Mi  ]  (  m3  ^ 

Rm  «•  \0,c,d)  =  l  J]  «".(<'  =  Repair) 1 |  (*’ :  No  "P3*1)  j 


(3.36) 


where  M\  =  number  of  members  being  repaired  or  replaced  at  the  j-th  inspection, 
M2  =  number  of  members  found  intact  at  the  j—\h  inspection,  and  M\  +  M2  =  M. 
In  Eq.  3.36,  Rm  (<*  :  Repair)  and  Rm  (t*  :  No  Repair)  are  identical  with  the  reliabilities 
R(t*  :  Repair)  and  R(t*  :  No  Repair)  defined  in  Eqs.  3.24  and  3.25,  respectively.  The  sub¬ 
script  m  is  used  to  indicate  that  these  reliabilities  are  associated  with  member  m.  Furthermore, 
the  dependence  of  Rm  (<*  :  Repair)  and  Rm  (t*  :  No  Repair)  as  implied  in  Eqs.  3.24  and  3.26 
has  not  been  explicitly  shown  for  the  sake  of  brevity. 

3.5.5  Time  Tj+i  for  (j+l)-th  Inspection 

If  the  reliability  of  the  entire  structure  is  specified  to  be  not  less  than  a  value  .^design’  ** 
such  that: 


Rm  (O  ^  ^design 


(3.37) 


Then  the  time  T}- ft  of  the  (j  4- 1)  —  th  inspection  is  found  as  the  minimum  value  of  t*  which 
satisfies  the  above  equation. 


3.6  Numerical  Example 

Numerical  simulations  have  been  carried  out  in  order  to  verify  the  validity  and  effectiveness 
of  Bayesian  analysis  to  determine  appropriate  inspection  intervals  and  true  values  of  uncertain 
parameters.  The  structure  used  in  this  study  is  assumed  to  have  100  members  (M  =  100).  Its 
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design  service  life  is  25  years  and  the  desired  minimum  reliability  level  of  the  entire  structure 
throughout  its  service  life  is  set  equal  to  0.80  (Resign  =  0.80).  In  general,  three  uncertain 
parameters  have  been  considered:  /?,  d  and  c.  Four  cases  have  been  worked  out  according  to 
the  number  of  uncertain  parameters  involved  in  each  one:  three  cases  involving  two  uncertain 
parameters  each  (Case  1:  (3  and  d  are  uncertain,  Case  2:  (3  and  c  are  uncertain  and  Case  3: 
d  and  c  are  uncertain)  and  one  case  involving  all  three  uncertain  parameters  (Case  4:  (3,  d  and  c 
are  uncertain).  The  true  values  of  the  uncertain  parameters  as  well  as  their  assumed  ranges 
are  shown  in  Table  3.1,  along  with  the  values  of  the  deterministic  parameters  appearing  in  the 
problem.  According  to  these  values,  the  probability  density  function  of  the  TTCI  is  plotted  in  Fig. 
3.3,  the  law  of  crack  propagation  in  Fig.  3.4,  the  probability  of  crack  detection  in  Fig.  3.5  and 
the  failure  rate  and  reliability  curves  before  and  after  craetc  initiation  in  Fig.  3.6.  As  an  example 
of  POD  curves  obtained  from  actual  data.  Fig.  3.7  displays  the  50%  (upper  curve)  and  95% 
(lower  curve)  confidence  level  curves  fitted  to  data  obtained  from  magnetic  particle  examination 
under  water  for  detection  of  fatigue  cracks  and  artificial  defects  (compare  to  Fig.  3.5).  Figure 
3.8  shows  the  failure  rate  after  crack  initiation  according  to  Ref.  3.2  (compare  to  Fig.  3.6). 

Before  dealing  with  the  above-mentioned  four  cases  involving  uncertain  parameters,  the 
particular  case  where  all  three  uncertain  parameters  (/ 3 ,  d  and  c)  assume  their  true  values  is 
examined.  The  inspection  schedule  resulting  from  one  simulation  is  displayed  in  Table  3.2 
while  the  corresponding  structural  reliability  for  the  entire  structure  is  plotted  in  Fig.  3.9.  The 
results  shown  in  Table  3.2  and  Fig.  3.9  are  obtained  in  the  following  way.  Equations  3.24 
and  3.25  are  used  to  calculate  Rm  (<*  |  (3,c,d )  given  by  Eq.  3.36.  In  this  case,  since  all  three 
uncertain  parameters  assume  their  true  values,  no  Bayesian  upgrading  is  needed.  Then,  the  time 
Ti+ 1  for  the  ( j  +  1)  -th  inspection  is  calculated  by  replacing  Rm  (t*)  appearing  in  Eq.  3.37 
with  Rm  ( t *  |  0,c,d).  Parameters  /3,c  and  d  appearing  in  the  expression  for  Rnf(t*  |  P,c,d) 
are  set  equal  to  their  true  values  shown  in  Table  3.1. 

Fifty  simulations  are  performed  for  each  of  the  four  cases  involving  uncertain  parameters. 
The  results  of  one  simulation  for  each  case  are  displayed  in  Tables  3.3-3.6.  These  results  include 
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the  inspection  schedule,  number  of  failed  members  and  number  and  length  of  detected  cracks. 
The  corresponding  structural  reliabilities  for  the  entire  structure  are  plotted  in  Figs.  3.10-3.13. 
The  results  shown  in  Tables  3.3-3.6  and  Figs.  3.10-3.13  are  obtained  in  the  following  way. 
Equations  3.24  and  3.2S  are  used  again  to  calculate  I  /?>c,  d)  given  by  Eq.  3.36.  In  this 

case,  since  uncertain  parameters  are  considered,  Bayesian  upgrading  is  needed.  This  upgrading 
is  carried  out  by  using  Eq.  3.34  to  compute  /J  (/?,  c,  d)  in  Eq.  3.35.  After  having  calculated 
Rm  (t*)  according  to  Eq.  3.35,  the  time  TJ+i  for  the  (j  +  1)  -th  inspection  is  calculated  using  Eq. 
3.37.  Finally,  the  corresponding  posterior  joint  density  functions  of  the  uncertain  parameters  after 
the  fifth  and  tenth  inspections  are  plotted  in  Figs.  3.14-3.16  for  cases  1,  2  and  3  respectively.  In 
these  figures,  it  can  be  easily  seen  that  for  the  specific  simulations  considered,  the  modal  values 
of  the  posterior  joint  density  functions  coincide  with  the  true  values  of  the  uncertain  parameters 
after  ten  inspections.  Unfortunately,  such  a  plot  cannot  be  given  for  case  4,  since  it  involves 
three  uncertain  parameters. 

The  statistics  of  the  modal  values  of  the  posterior  joint  density  functions  after  five  and  ten 
inspections  are  shown  in  Table  3.7  for  the  four  cases  examined.  Finally,  the  average  inspection 
frequency  during  the  twenty-five  years  of  design  service  life  is  displayed  in  Table  3.8  for  the 
four  cases  examined.  It  should  be  noted  that  the  statistics  appearing  in  Tables  3.7  and  3.8  are 
calculated  from  the  results  of  the  fifty  simulations  performed. 


3.7  Structures  With  Members  Subjected  To  Different  Stress  Levels 
3.7.1  Introduction 

The  structures  considered  in  Sections  3.2-3.6  consisted  of  structural  components  subjected 
to  the  same  level  of  stress  intensity.  In  this  chapter,  however,  the  structures  are  assumed  to 
consist  of  three  classes  of  components:  class  A  components  subjected  to  the  highest  stress  level 
Ay  class  B  components  subjected  to  stress  level  B  and  class  C  components  subjected  to  stress 
level  C  where  A  >  B  >  C. 
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3.7.2  Bayesian  Analysis 


The  reliability  of  the  group  of  structural  components  subjected  to  stress  intensity  level  J  at 


38) 


R 


time  i*  after  the  j— th  inspection  is  denoted  by  Rm,j  (<*)  and  is  given  by: 

fijma,  dmas  CJmas 

757(0-  J  J  J  (3. 

ftjmin  dmin  CJmin 

where  J  =  A,B,C  and: 

{Mu  ]  (Mu  } 

n  «mj(C:Repair)|  |  J]  RmJ  (f:No  Repair)  j  (3.39) 

where  M\j  =  number  of  members  of  group  J  being  replaced  or  repaired  at  the  j— th  inspection, 

M2j  =  number  of  members  of  group  J  found  intact  at  the  j -th  inspection,  and  M\j+M2j  =  Mj 

=  number  of  members  in  group  J.  The  posterior  joint  density  function  of  the  three  uncertain 

parameters  of  group  J  immediately  after  the  j— th  inspection  is  given  by: 

r  j  /  a  j  \ _ LFjj  x  fj  cj) _ 

Ji  i  J)  j  j  j  (Numerator)  dfljd  (d)  dcj 
Dj 

Therefore,  the  structural  reliability  of  the  entire  structure  can  be  estimated  as: 


(3.40) 


Ru(n=  n 

J=A,B,C 

Time  r;+1  of  the  ( j  +  l)-th  inspection  can  then  be  estimated  with  the  aid  of  Eq.  3.37. 


(3.41) 


3.73  Parameter  Values 


The  numerical  values  of  the  various  parameters  appearing  in  the  analytical  model  are  assumed 
to  be  the  same  as  those  chosen  in  the  previous  chapter,  with  the  obvious  exception  of  parameters 
that  depend  on  the  stress  intensity  level. 

In  the  numerical  analysis  that  follows,  it  is  assumed  that  the  stress  intensity  characteristics 
of  the  three  groups  of  structural  components  are  given  by: 


Stress  intensity  level 

A 

B 

c 

Stress 

S 

0.9S 

0.8S 

Stress  range 

Act 

0.9A<t 

0.8A<7 
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On  the  basis  of  the  above  table,  the  following  assumptions  are  made  for  those  parameter 
values  that  depend  on  the  stress  intensity  level: 

a.  Scale  Parameter  P  in  TTCI  Distribution: 

If  P  =  30  years  under  a  cyclic  load  with  amplitude  Act,  the  values  of  P  for  cyclic  loads  with 
amplitudes  0.9A a  and  0.8Acr  are  estimated  to  be  45  years  and  65  years  respectively: 


Stress  intensity  level 

A 

B 

C 

P  (years) 

30 

45 

65 

b.  Coefficient  c  in  Crack  Propagation  Law: 

According  to  Eq.  3.3, 

~  =  Ci  (A K)b  =  Ci  (Ac 7y/^)b  =  c  •  a  (3.42) 

where  6  =  2  and  c  =  7rCi  (Act)2  .  Therefore,  c  is  proportional  to  (Act)2  and  hence, 


Stress  intensity  level 

A 

B 

c 

c 

0.6 

issBEna 

c.  Parameters  r  and  q  in  Failure  Rate  Expressions: 

To  be  rigorous,  the  values  of  parameters  r  and  q  for  different  stress  levels  must  be  obtained 
by  evaluating  the  crossing  rate  of  appropriate  threshold  values  of  the  response  stress  process. 
Since  in  the  present  study  no  specific  dynamic  system  is  considered  for  a  response  analysis 
and  since  emphasis  is  placed  on  the  development  of  a  Bayesian  method  for  nondestructive 
inspection  procedures,  the  values  shown  in  the  following  table  are  assumed  for  r  and  q. 
Note  that  r  is  proportional  to  the  stress  intensity  level  and  q  is  proportional  to  its  square. 


65 


Stress  intensity  level 

A 

B 

c 

r 

-7.5 

-9.0  =  1.2  x  (-7.5) 

q 

0.9 

H22JEEGEEEI 

0.576  =  0.82  x  0.9 

The  parameter  values  for  0,  c,  r,  and  q  shown  in  the  previous  three  tables  will  be  used  for 
Monte  Carlo  simulation  study. 

3.7.4  Numerical  Example 

Numerical  simulations  have  been  carried  out  in  order  to  verify  the  validity  and  effectiveness 
of  Bayesian  analysis  to  determine  inspection  intervals  and  true  values  of  uncertain  parameters. 
The  structure  used  in  this  study  is  assumed  to  have  100  members  (Mr  -  100).  Of  the  total 
number  of  100  members,  20  are  subjected  to  stress  intensity  level  A,  30  to  stress  intensity  level 
B  and  50  to  stress  intensity  level  C. 

The  design  service  life  is  25  years  and  the  desired  minimum  reliability  level  of  the  entire 
structure  throughout  its  service  life  is  set  equal  to  0.80  (.RdeSjgn  =  0.80). 

Only  one  case  involving  two  uncertain  parameters  (specifically  0  and  c)  is  examined  now. 
The  true  values  of  the  uncertain  parameters  as  well  as  their  assumed  ranges  are  shown  in  Table 
3.9,  along  with  the  values  of  the  deterministic  parameters  appearing  in  the  problem. 

Before  dealing  with  the  above-mentioned  case  involving  two  uncertain  parameters,  the 
particular  case  where  both  uncertain  parameters  (0  and  c)  assume  their  true  values  is  examined. 
The  inspection  schedule  resulting  from  one  simulation  is  displayed  in  Table  3.10,  while  the 
corresponding  structural  reliability  for  the  entire  structure  and  for  the  three  stress  intensity  level 
groups  is  plotted  in  Fig.  3.17. 

Fifty  simulations  are  now  performed  for  the  case  involving  the  two  uncertain  parameters 
0  and  c.  The  results  of  one  of  these  simulations  are  displayed  in  Table  3.11.  These  results 
include  the  inspection  schedule,  number  of  failed  members  in  each  stress  intensity  level  group, 
and  number  and  length  of  detected  cracks  in  each  stress  intensity  level  group.  The  corresponding 
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structural  reliabilities  for  the  entire  structure  and  for  the  three  stress  intensity  level  groups  are 
plotted  in  Fig.  3.18.  Finally,  the  corresponding  posterior  joint  density  functions  of  the  two 
uncertain  parameters  after  the  third  and  sixth  inspection  are  plotted  in  Fig.  3.19  for  the  three  stress 
intensity  level  groups,  A,  B ,  and  C.  The  estimation  of  the  true  values  of  uncertain  parameters  /? 
and  c  immediately  after  the  sixth  inspection  for  structural  components  under  stress  intensity  level 
A  is  accomplished  reasonably  well.  This  is  due  primarily  to  the  fact  that  8  cracks  were  found 
during  the  first  six  inspections.  It  should  be  noted  that  die  estimation  is  considered  reasonable 
in  the  sense  that  a  distinct  mode  emerges  at  a  location  close  to  the  true  values  of  /?  and  c.  On 
the  other  hand,  the  same  did  not  apply  to  structural  components  under  stress  intensity  levels 
B  and  C  for  which  only  7  and  3  cracks  were  respectively  found.  In  all  cases,  however,  the 
posterior  probability  density  function  after  the  sixth  inspection  indicates  a  mode  located  closer 
to  the  true  location  than  after  the  third  inspection. 

The  statistics  of  the  modal  values  of  the  posterior  joint  density  functions  after  three,  six  and 
seven  inspections  are  shown  in  Table  3.12.  These  statistics  are  based  on  Monte  Carlo  simulations 
with  a  sample  size  equal  to  SO.  Finally,  the  average  number  of  inspections  required  to  maintain 
the  specified  reliability  level  for  25  years  is  equal  to  6.8,  based  on  Monte  Carlo  simulations  with 
a  sample  size  equal  to  50.  This  value  of  6.8  should  be  compared  with  6.0  which  is  die  number 
of  inspections  required  when  the  true  values  of  /?  and  c  are  known  (see  Table  3.13). 

3.8  Future  Work 

The  following  five  areas  require  further  study: 

a.  In  this  work,  the  three  uncertain  parameters  /?,c,  and  d  were  considered  to  be  uncorrc- 
iated.  However,  there  is  strong  evidence  that  /?,  c,  and  d  are  in  reality  correlated  to  each 
other.  Therefore,  the  statistical  correlation  among  0,c,  and  d  and  the  influence  of  this 
correlation  on  the  obtained  results  require  further  study.  Another  part  of  future  work 
concerning  parameters  0,  c,  ami  d  is  to  examine  their  sensitivity  to  the  obtained  results; 
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b.  The  effect  of  the  form  of  certain  POD  curves  on  the  reliability  of  marine  structures 
throughout  their  service  life  requires  further  study.  A  comparison  has  to  be  made  among 
several  established  POD  curves  in  the  aerospace  industry,  in  order  to  assess  their  relative 
influence  on  the  reliability  of  marine  structures  subjected  to  non-periodic  inspections.  In 
this  way,  more  reliable  POD  curves  can  be  established  for  marine  structures; 

c.  The  cost-effectiveness  of  the  proposed  method  of  non-periodic  inspections  based  on 
Bayesian  analysis  requires  further  study.  Specifically,  a  cost-benefit  analysis  can  be 
performed  taking  into  consideration  the  cost  of  the  non-periodic  inspection  procedure 
and  the  increased  level  of  reliability  for  the  structure.  These  results  have  to  be  compared 
with  the  results  of  cost-benefit  analysis  associated  with  die  standard  periodic  inspection 
procedure; 

d.  Verification  of  the  proposed  methodology  using  actual  data  from  inspections  of  marine 
structures  is  one  of  the  most  important  tasks  of  future  work.  This  task  can  be  accom¬ 
plished  by  taking  advantage  of  already  completed  inspections  in  marine  structures  to 
determine  whether  these  structures  actually  maintained  a  prespecified  reliability  level 
throughout  their  life; 

e.  The  failure  rate  expression  after  crack  initiation  should  at  least  be  validated  by  Monte 
Carlo  simulation  utilizing  the  crack  propagation  law  and  uncontrolled  crack  growth 
condition  based  on  fracture  mechanics  theory  under  various  random  stress  histories 
consistent  with  the  stress  intensity  factor  fluctuation  considered. 
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Table  3.1  Parameter  Values  of  Numerical  Example, 


Item 

True  value 

Assumed 

range 

Design  life 

25  years 

— 

Minimum  structural  reliability 
required 

R,  < 
design 

0.8 

Total  number  of 
structural  members 

M 

100 

— — 

Parameters  in  PDF  of  TTCI 

o 

1  0 

2.0 

30  years 

20-40 

Parameter  in  crack 
propagation 

c 

0.6 

0.4  -0.6 

Initial  crack  length 

a0 

10  mm  (  0.4  in  ) 

- — 

Parameter  in  POD 

d 

0.01 

0.002-0.018 

Parameters  in  failure  rate 

<5 

-7.5 

0.9 

— 

Table  3.2 


Inspection  Schedule  for  the  Case  Where  All  Uncertain  Parameters  Assume  Their 
True  Values  (0  =  30,  d  =  0.01  and  e  =  0.6). 


Inspection 

No. 

Inspection 
time  :  T 
(years) 

1 

4.9 

2 

7.9 

3 

10.1 

4 

12.1 

5 

13.9 

6 

15.7 

7 

17.4 

8 

19.0 

9 

20.6 

10 

22.1 

11 

23.6 
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Table  3.3  Results  for  Case  1  (Uncertain  Parameters  :  0  and  d) 


Xnapectlon 

No. 

Inspection 
tine  :  T 
(years) 

No.  of 
failed 
aenbars 

No.  of 

detected 

cracks 

Detected  crack  length  (cm) 

1 

S.O 

1(b) 

1 

22 

2 

7.9 

0 

1 

42 

3 

10.1 

0 

4 

189,189,124,  23 

4 

12.1 

0 

1 

1215, 

5 

13.5 

0 

3 

117,  29,  25 

6 

15.2 

0 

4 

110,110,  54,  98 

7 

17.0 

0 

7 

16,241,  29,  13,465.  98,  21 

8 

18.7 

1(a) 

6 

64,178,168.104,  42,117 

9 

20.3 

0 

4 

68,124,241,  45 

10 

11.9 

0 

3 

149,140,124 

11 

23.3 

0 

4 

628,189,227,158 

12 

24.9 

0 

6 

168,  98,  64,  57,  42,306 

(a)  :  Failure  after  crack  Initiation  1  in  -  25.4  mm 

(b)  t  Failure  before  crack  Initiation 


Table  3.4  Results  for  Case  2  (Uncertain  Parameters  :  0  and  e) 


Inspection 

No. 

Inspection 

time  :  T 
(years) 

No.  of 
failed 
members 

No.  of 

detected 

cracks 

! 

Detected  crack  length  (mm) 

1 

5.0 

1(b) 

1 

22 

2 

l  7.8 

0 

1 

40 

3 

9.6 

0 

3 

140,140.  92 

4 

11.2 

0 

1 

708 

5 

13.4 

0 

3 

33,168,  11 

6 

15.3 

0 

4 

117,345,  87,  87 

7 

17.2 

0 

8 

366.178,  72,  33,  14,227,525,  23 

8 

19.0 

1(a) 

7 

42,201,124,110,  51,  21,140 

9 

20.5 

0 

5 

189,140.  26,  61,  54 

10 

22.0 

0 

3 

189,  29,  54 

11 

23.5 

0 

4 

798,189,227,389 

(a)  :  Failure  after  crack  Initiation  1  in »  25.4  mm 

(b)  :  Failure  before  crack  initiation 
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Table  3.5  Results  for  Case  3  (Uncertain  Parameters  :  d  and  c). 


Inspection 

No. 

Inspection 
time  :  T 
(years) 

m 

No.  of 

detected 

cracks 

Detected  crack  length  (mo) 

1 

4.9 

1(b) 

1 

21 

2 

7.9 

0 

1 

42 

3 

9.8 

0 

3 

158,158,104 

4 

11.6 

0 

1 

900 

5 

13.2 

0 

3 

21,149,  48 

6 

15.0 

0 

4 

288,  98,  48,  87 

7 

16.8 

0 

5 

288,213,213,  26,  87 

8 

18.6 

Ka) 

8 

61,168,158,  98,  87,  40,  16,110 

9 

20.4 

0 

5 

98,132,  98,  57,158 

10 

22.2 

0 

6 

213,366,  87,178,  72,158 

11 

23.9 

0 

3 

54,288,168 

(a)  :  Failure  after  crack  initiation  1  in  -25.4  mm 

(b)  :  Failure  before  crack  initiation 
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Table  3.6  Results  for  Case  4  (Uncertain  Parameters  :  /?,  d  and  c). 


Inspection 

No. 

Inspection 
time  :  T 
(years) 

No.  of 
failed 
members 

No.  of 

detected 

cracks 

Detected  crack  length  (mm) 

1 

5.0 

1(b) 

1 

22 

2 

7.8 

0 

1 

40 

3 

9.6 

0 

3 

140,140,  92 

4 

11.2 

0 

1 

708 

5 

12.9 

0 

1 

124 

6 

14.7 

0 

5 

82,241,  82,  61,  61 

7 

16.6 

0 

8 

124,  51,227,158,366,  77,  12,  33 

8 

18.4 

1(a) 

4 

37,140,  87,  48 

9 

20.1 

0 

6 

82,149,110,213,  82,  42 

10 

21.8 

0 

6 

168,  26,  68,  82,525,110 

11 

23.5 

0 

5 

42,389,158,132,366 

(a)  :  Failure  after  crack  Initiation  lin»25.4mra 

(b)  :  Failure  before  crack  initiation 
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Table  3.7  Statistics  of  Modal  Values  of  the  Posterior  Joint  Density  Functions  of  Uncertain 
Parameters  (True  Values  :  0  —  30,  d  =  0.01  and  c  =  0.6;  Number  of  Simulations 
=  50) 


1.  The  Fifth  Inspection 


Mean 


Case 

6 

d 

c 

1 

31.6 

0.0112 

— 

2 

32.1 

— 

0.622 

3 

— 

0.0109 

0.616 

4 

32.2 

0.0121 

0.617 

Standard 

deviation 

1 

Case 

6 

d 

c 

1 

4.23 

0.000354 

— 

2 

4.10 

— 

0.0884 

3 

— 

0.000361 

0.0924 

4 

4.47 

0.000382 

0.0976 

C0V 


Case 

8 

d 

c 

1 

0.134 

0.292 

— 

2 

0.128 

— 

0.142 

3 

— 

0.331 

0.150 

4 

0.139 

0.316 

0.158 

2.  The  Tenth  Inspection 


Mean 


Case 

6 

d 

c 

1 

31.4 

0.0114 

— 

2 

31.1 

— 

0.606 

3 

— 

0.0116 

0.635 

4 

31.5 

0.0122 

0.636 

Standard 

deviation 

Case 

8 

d 

c 

1 

2.83 

0.000238 

— 

2 

2.84 

— 

0.0646 

3 

— 

0. 000280 

0.0801 

4 

2.92 

0.000303 

0.0845 

COV 


Case 

8 

d 

c 

i 

0.0902 

0.209 

— 

2 

0.0913 

■ - 

0.107 

3 

— 

0.241 

0.126 

4 

0.0926 

0.248 

0.133 
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Table  3.8 


Average  Inspection  Frequency  During  25  years  (Number  of  Simulations  =  50) 


Case 

Frequency 

1 

11.4 

2 

11.6 

3 

11.6 

4 

11.4 

True 

value 

11 
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Item 

Model  value 

Design  life  (years) 

25 

Minimum  structural  reliability  design 
Required 

0.8 

Stress  intensity  level 

A 

B 

C 

Stress 

S 

0.9S 

0.8S 

Stress  range 

6a 

0.9  Ac 

0.8  6a 

Numl.er  of  structural  members  M 
(Total  :  Mt-100) 

20 

30 

50 

Parameters  in  PDFof  TTCI 

a 

2.0 

2.0 

2.0 

0  (years)  True  value 

Assumed  range 

30 

(20-40) 

45 

(35-55) 

65 

(55-75) 

Parameter  in  crack  propagation 

b 

2.0 

2.0 

2.0 

c  True  value 

Assumed  range 

0.8 

(0.4 -0.8) 

0.486 

(0.286-0.886) 

0.384 

(0.184-0.584) 

Initial  cracklength  aQ  mm(in) 

10(0.4) 

10(0.4) 

10(0.4) 

Parameter  in  POD  d 

0.01 

0.01 

0.01 

Parameter  in  failure  rate  r 

-7.5 

•8.25 

-9 

q 

0.9 

0.729 

0.576 

Table  3.9  Parameter  Values  of  Numerical  Example 

(  Case  2  Uncertain  Parameters  :  0  and  c  ] 
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Inspection 

No. 

Inspection 
time  :  T 
(years) 

1 

8.3 

2 

12.1 

3 

15.3 

4 

18.2 

5 

20.9 

6 

23.6 

Table  3.10  Inspection  Schedule  for  True  Value 


Inspection  »«•  ' »*  Detected  crack  length 

time  :  T  failed  detected  ,  »  & 

l,°*  (years)  members  cracks  v  ; 
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The  Third  Inspection  2.  The  Sixth  Inspection  3.  The  Seventh  Inspection 


Table  3.12  Statistics  of  Modal  Value  for  Posterior  Joint  Probability  Density  of 
(  Case  2  Uncertain  Parameters  :  0  and  c  ) 

f  A(30,  0.6) 

True  Value  (0,  c)  -  B(45,  0.486) 

1  C(65,  0.384) 


-fr,:  •'  ~  ■:■ 


Case 

Frequency 

2 

6.8 

True 

value 

6 

Table  3.13  Average  Inspection  Frequency 
during  25  years 

,  Case  2  Uncertain  Parameters  : 
1  Number  of  Simulations  *  50 


6  and  c  j 


CRACK  LENGTH 


Fig.  3.1  Total  Fatigue  Life  Consists  of  Crack  Initiation,  Subcriti- 
cal  Crack  Propagation  and  Final  Fracture. 
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Ti 
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crack  initiates 
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Fig.  3.2b 


Ti 
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Ti+i 
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i- 1 


interval  where 
failure  occurs 


Fig.  3.2c 
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*V-t 
H — 


crack  initiates 
here 


interval  where 
failure  occurs 


Fig.  3.2d 


Ti 
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Ti 
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tc 


Ti- 1 

H — 


r< 


crack  initiates 
here 


t  Ti+ 1 

■i — I — 


*>-i 

H — 


crack  initiates 
here 


Fig.  3.2e 


Fig.  3.2f 


Fig.  3.2  Chronology  of  Events. 
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Fig.  3.3  Probability  Density  Function  of  Time  to  Crack  Initiation. 
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Crack  length  :  a(t)  (mm) 


vjpip,4»!P 


ifUMPiiimsiMi 


Fig.  3.4  Law  of  Crack  Propagation. 


86 


I  100 

a0  *  10 


300  400  500 

Crack  length  :  a  (mm) 


Fig.  3.5  Probability  of  Crack  Detection. 
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Fig.  3.7  50%  (upper  curve)  and  95%  (lower  curve)  Confidence 
Level  Curves  Fitted  to  Data  Obtained  from  Magnetic 
Particle  Examination  Applied  Under  Water  for  Detection 
of  Fatigue  Cracks  and  Artificial  Defects. 


FAILURE  RATE  HIT1 


Fig.  3.8  Failure  Rate  After  Crack  Initiation  According  to  [3.2]. 
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0 


J  5 


10 


15  20  25 

Service  time  :  t  (years) 


Fig.  3.9  Inspection  Schedule  and  Structural  Reliability  for  the 
Case  Where  All  Uncertain  Parameters  Assume  Their 
True  Values. 


0 


5 


10 


15  20  25 

Service  time  :  t  (years) 


Fig.  3.10  Inspection  Schedule  and  Structural  Reliability  for  Case  1 
(Uncertain  Parameters:  (3  and  d). 
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Reliability  Rm 


True 


Service  time  :  t  (years) 


Fig.  3.11  Inspection  Schedule  and  Structural  Reliability  for  Case  2 
(Uncertain  Parameters:  (3  and  c). 
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Reliability 


True 


Fig.  3.12  Inspection  Schedule  and  Structural  Reliability  for  Case  3 
(bncertain  Parameters:  d  and  c). 
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Reliability  Rm 


True 


Service  time  :  t  (years) 


Fig.  3.13  Inspection  Schedule  and  Structural  Reliability  for  Case  4 
(Uncertain  Parameters:  (3 ,  d  and  c). 
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Posterior  probability  Posterior  probability 

density  f(0,d)  density  f($,d) 


20 


1.  The  Fifth  Inspection 


True  value  (3,  d)  -  (30,  0.01) 


Fig.  3.14  Posterior  Joint  Probability  Density  Function  for  Case  1 
(Uncertain  Parameters:  and  d). 
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Posterior  probability 


1.  The  Fifth  Inspection 


True  value  (B,  c)  ■  (30,  0.6) 


Fig.  3.15  Posterior  Joint  Probability  Density  Function  for  Case  2 
(Uncertain  Parameters:  /?  and  c). 


Posterior  probability  Posterior  probability 

density  f(c,d)  density  f(c,d) 


0.5 


(cM,  dM)-  (0.6,  0.00733) 


1.  The  Fifth  Inspection 


True  value  (c,  d)  ■  (0.6,  0.01) 


Fig.  3.16  Posterior  Joint  Probability  Density  Function  for  Case  3 
(Uncertain  Parameters:  d  and  c). 
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Rellabilit 


0  5  10  15  20  25 

Service  time  :  t  (years) 


Fig.  3.17a  Inspection  Schedule  and  Structural  Reliability  for  True 
Value  (Whole  Structure). 
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Reliability  R^ 


Fig.  3.17b  Inspection  Schedule  and  Structural  Reliability  for  True 
Value  (Stress  Intensity  Level  A). 
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Reliability  R^r 


I  ....  i  ■ _ i,l  ■ — i — . — i — L 

0  5  10  15  20  25 

Service  tiuie  :  t  (years) 


Fig.  3.17d  Inspection  Schedule  and  Structural  Reliability  for  True 
Value  (Stress  Intensity  Level  C). 
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10 


■  -  - _ . _ . _ I _ I—J - 1 - 1 - L. 

15  20  25 

Service  time  :  t  (years) 


Fig.  3.18a  Inspection  Schedule  and  Structural  Reliability  for  Case 
2:  Uncertain  Parameters  f3  and  c  (Whole  Structure). 
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Reliability 


0  5  10  15  20  25 

Service  time  :  t  (years) 


Fig.  3.18b  Inspection  Schedule  and  Structural  Reliability  for  Case 
2:  Uncertain  Parameters  and  c  (Stress  Intensity  Level 
A). 
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Reliability  Rmb 


True 


Fig.  3.18c  Inspection  Schedule  and  Structural  Reliability  for  Case 
2:  Uncertain  Parameters  (3  and  c  (Stress  Intensity  Level 

B). 
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Reliability  R^c 


True 


0.8 


L  ■  ■  -  i  l  1—  . - 1  -  i  l  I  i  i  «  ,<  X- 

0  5  10  15  20  25 

Service  tine  :  t  (years) 


Fig.  3.18d  Inspection  Schedule  and  Structural  Reliability  for  Case 
2:  Uncertain  Parameters  0  and  e  (Stress  Intensity  Level 
C). 
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1.  The  Third  Inspection 


2.  The  Sixth  Inspection 


True  value  (B,  c)  *  (30,  0.6) 


Fig.  3.19a  Posterior  Joint  Probability  Density  Function  for  Case  2: 

Uncertain  Parameters  0  and  c  (Stress  Intensity  Level  A). 
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0.6 


-(6m,  cM)  -  (35,  0.419) 


1.  The  Third  Inspection 


•H  - 
J3  V 
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U  to 
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in  *a 
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a. 


2.  The  Sixth  Inspection 


True  value  (6,  c)  ■  (45,  0.486) 


Fig.  3.19b  Posterior  Joint  Probability  Density  Function  for  Case  2: 

Uncertain  Parameters  and  c  (Stress  Intensity  Level  B). 
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1.  The  Third  Inspection 


0.584, 


2.  The  Sixth  Inspection 

True  value  (8,  c)  ■  (65,  0.384) 


Fig.  3.19c  Posterior  Joint  Probability  Density  Function  for  Case  2: 

Uncertain  Parameters  /?  and  c  (Stress  Intensity  Level  C). 
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